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Abstract In this lecture the basic theory of accretion disks is reviewed, with emphasis on 
aspects relevant for X-ray binaries and Cataclysmic Variables. The text gives a general intro- 
duction as well as a selective discussion of a number of more recent topics. 



p:^ 1 Introduction 

(— I Accretion disks are inferred to exist as objects of very different scales: millions of km in 

I low Mass X-ray Binaries (LMXB) and Cataclysmic Variables (CV), solar radius-to-AU scale 

JL disks in protostellar objects, and AU-to-parsec scale disks in Active Galactic Nuclei (AGN). 
%^ An interesting observational connection exists between accretion disks and jets (such as 

c/3 the spectacular jets from AGN and protostars), and outflows (the 'CO-outflows' from proto- 

Cd stars and the 'broad-line-regions' in AGN). Lacking direct (i.e. spatially resolved) observa- 
tions of disks, theory has tried to provide models, with varying degrees of success. Uncer- 

1—1 tainty still exists with respect to some basic questions. In this situation, progress made by 

^ observations or modeling of a particular class of objects has direct impact for the understand- 

^^ ing of other objects, including the enigmatic connection with jets. 

Lj In this lecture I concentrate on the more basic aspects of accretion disks, but an attempt 

l^ is made to mention topics of current interest as well. Some emphasis is on those aspects of 

• accretion disk theory that connect to the observations of LMXB and CV's. For other reviews 

J— ^ on the basics of accretion disks, see Pringle (1981), Papaloizou and Lin (1995). For a more 

f^ extensive introduction, the textbook by Frank et al. (2002). For a comprehensive text on CVs, 

1— I Warner 1995. 

> 

X 

%-t 2 Accretion: general 

Gas falling into a point mass potential 

GM 

= 

r 
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from a distance ro to a distance r converts gravitational into kinetic energy, by an amount 
A0 = GM{\/r— 1/ro). For simplicity, assuming that the starting distance is large, A4> = 
GM/r. The speed of arrival, the free-fall speed Vff is given by 

^vl = GM/r. (1) 

If the gas is then brought to rest, for example at the surface of a star, the amount of energy e 
dissipated per unit mass is 

I 2 GM 

If, instead, it goes into a circular Kepler orbit at distance r: 

IGM , ,., 

e = [orbit ). 

2 r 

The dissipated energy may go into internal energy of the gas, and into radiation which escapes 
to infinity (usually in the form of photons, but neutrino losses can also play a role in some 
cases). 



2.1 Adiabatic accretion 

Consider first the case when radiation losses are neglected. Any mechanical energy dissipated 
stays locally in the flow. This is called an adiabatic flow (not to be confused with isentropic 
flow). For an ideal gas with constant ratio of specific heats 7, the internal energy per unit mass 
is 



(7-l)P 
With the equation of state 

P^MpT/ix (2) 

where ^ is the gas constant and fi the mean atomic weight per particle, we find the tempera- 
ture of the gas after the dissipation has taken place (assuming that the gas goes into a circular 
orbit): 

T=\{Y-l)Tv.. (3) 

where Tyir, the virial temperature is defined as 

where g is the acceleration of gravity at distance r. In an atmosphere with temperature near 
Tyir, the sound speed c^ — {jMT / nY''^ is close to the escape speed from the system, and the 
hydrostatic pressure scale height, H = ^T /{ng) is of the order of r. Such an atmosphere may 
evaporate on a relatively short time scale in the form of a stellar wind. This is as expected 
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from energy conservation: if no energy is lost through radiation, the energy gained by the 
fluid while falling into a gravitational potential is also sufficient to move it back out again. 

A simple example is spherically symmetric adiabatic accretion (Bondi, 1952). An im- 
portant result is that such accretion is possible only if 7 < 5/3. The lower 7, the lower the 
temperature in the accreted gas (eq. l3]l, and the easier it is for the gas to stay bound in the 
potential. A classical situation where adiabatic and roughly spherical accretion takes place is 
a supernova implosion: when the central temperature becomes high enough for the radiation 
field to start disintegrating nuclei, 7 drops and the envelope collapses onto the forming neu- 
tron star via an accretion shock. Another case are Thorne-Zytkow objects (e.g. Cannon et al. 
1992), where 7 can drop to low values due to pair creation, initiating an adiabatic accretion 
onto the black hole. 

Adiabatic spherical accretion is fast, taking place on the dynamical time scale: something 
on the order of the free fall time scale, or Kepler orbital time scale, 

Td = r/vK = i2K= {r^/GMfl^, (4) 

where vk, X2k are the Kepler orbital velocity and angular frequency. 

When radiative loss becomes important, the accreting gas can stay cool irrespective of the 
value of 7, and Bondi's critical value 7=5/3 plays no role. With such losses, the temperatures 
of accretion disks are usually much lower than the virial temperature. 



2.2 Temperature near compact objects 

For accretion onto a neutron star surface, R = \Q km, M = lAM^, we have a free fall 
speed Vff/c « 0.4 c (this in Newtonian approximation, the correct value in General Rela- 
tivity is quantitatively somewhat different). The corresponding virial temperature would be 
Ty ^2 10'^ K, equivalent to an average energy of 150 MeV per particle. 

This is not the actual temperature we should expect, since other things happen before 
such temperatures are reached. If the accretion is adiabatic, one of these is the creation of 
a very dense radiation field. The energy liberated per infalling particle is still the same, but 
it gets shared with a large number of photons. At temperatures above the electron rest mass 
(« 0.5Mey), electron-positron pairs e* are produced in addition to photons. These can take 
up most of the accretion energy, limiting the temperature typically to a few MeV. 

In most observed disks temperatures do not get even close to an MeV, however, because 
accretion is rarely adiabatic. Energy loss takes place by escaping photons (or under more 
extreme conditions: neutrinos). Exceptions are the radiatively inefficient accretion flows dis- 
cussed in section[T3] 



2.2.1 Radiative loss 

Next to the adiabatic temperature estimates, a useful characteristic number is the 'black body 
effective temperature'. Here, the approximation made is that the accretion energy is radiated 
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away from an optically thick surface of some geometry, under the assumption of a balance 
between the heating rate by release of accretion energy and cooling by radiation. For a specific 
example consider the surface of a star of radius R and mass M accreting via a disk. Most of the 
gravitational energy is released close to the star, in a region with surface area of the order, let's 
call it 4nR^. If the surface radiates approximately as a black body of temperature T (make 
a note of the fact that this is a bad approximation if the opacity is dominated by electron 
scattering), the balance would be 

M^^4nR'a,T\ (5) 

R 

where Cr is the Stefan-Boltzmann constant, Cr = a^c/A if a^ is Planck's radiation constant. 
For a neutron star with M = IAMq, = 3 10^^ E, R = 10 km, accreting at a typical observed 
rate (near Eddington rate, see below), M = lO'^g/s w lO^^Mp/yr, this temperature would be 
T w 10^ K, or « 1 keV per particle. Radiation with this characteristic temperature is observed 
in accreting neutron stars and black holes in their so-called 'soft X-ray states' (as opposed to 
their 'hard' states, in which the spectrum is very far from a black body). 

The kinds of processes involved in radiation from accretion flows form a large subject 
in itself and will not be covered here. For introductions see Rybicki & Lightman (1979) 
and Frank et al. (2002). At the moderate temperatures encountered in protostars and white 
dwarf accreters, the dominant processes are the ones known from stellar physics: molecular 
and atomic transitions and Thomson scattering. Up to photon energies around 10 keV (the 
Lyman edge of an iron nucleus with one electron left) these processes also dominate the 
spectra of neutron star and black hole accreters. Above this energy, observed spectra become 
dominated by Compton scattering. Cyclotron/synchrotron radiation plays a role in when a 
strong magnetic field is present, which can happen in most classes of accreting objects. 



2.3 Critical luminosity 

Objects of high luminosity have a tendency to blow their atmospheres away due to the radia- 
tive force exerted when the outward traveling photons are scattered or absorbed. Consider a 
volume of gas on which a flux of photons is incident from one side. Per gram of matter, the 
gas presents a scattering (or absorbing) surface area of K cm^ to the escaping radiation. The 
force exerted by the radiative flux F on one gram is Fk/c. The force of gravity pulling back 
on this one gram of mass is GM/r^. The critical flux at which the two forces balance (energy 
per unit area and time) is 

Fe^ J-- (6) 

Assuming that this flux is spherically symmetric, it can be converted into a luminosity, 

LE = 47rGMc/K-, (7) 

the Eddington critical luminosity, popularly called the Eddington limit (e.g. Rybicki and 
Lightman, 1979). If the gas is fully ionized, its opacity is dominated by electron scatter- 
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ing, and for solar composition K is then of the order 0.3 cm^/g (about a factor 2 lower for 
fully ionized heUum). With these assumptions, 

This number is different if the opacity is not dominated by electron scattering. In partially 
ionized gases of solar composition bound-bound and bound-free transitions can increase the 
opacity by a factor up to lO"*; the Eddington flux is then correspondingly lower 

If the luminosity results from accretion, one can define a corresponding Eddington char- 
acteristic accretion rate Me: 

GM . . , 

Me = Le -^ ME = 47trc K. (8) 

r 

With K = 0.3: 

Me « 1.310^^6 g/s^llO-^re M^yr-^ 

where rg is the radius of the accreting object in units of 10*' cm. The characteristic accretion 
rate thus scales with the size of the accreting object, while the critical luminosity scales with 
mass. 

Whereas Le is a critical value which in several circumstances plays the role of a limit, the 
Eddington characteristic accretion rate is less of a limit. For more on exceptions to Le and 



Me see 2.3.2 below. 



2.3.1 Eddington luminosity at high optical depth 

The Eddington characteristic luminosity was derived above under the assumption of a radi- 
ation flux passing through an optically thin medium surrounding the radiation source. What 
changes if the radiation passes through an optically thick medium, such as a stellar interior? 
At high optical depth the radiation field can be assumed to be nearly isotropic, and the diffu- 
sion approximation applies (cf. Rybicki & Lightman 1979). The radiative heat flux can then 
be written in terms of the radiation pressure Pi- as 

"■"§' <9) 

dr 

where T is the optical depth 

dT = K-pd,s (10) 

along a path s, and K an appropriate frequency-averaged opacity (such as the Rosseland 
mean). Balancing the gradient of the radiation pressure against the force of gravity gives 
the maximum radiation pressure that can be supported: 

VPr,max=gP, (H) 
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where g is the acceleration of gravity. With (|9| this yields the maximum radiation flux at a 
given point in a static gravitating object: 

Fr.max = Cg/KT, (12) 

i.e. the same as the critical flux in the optically thin case. 

2.3.2 Limitations of the Eddington limit 

The derivation of Fe assumed that the force relevant in the argument is gravity. Other forces 
can be larger. An example would be a neutron star with a strong magnetic field. The curvature 
force B^/{4nrc) in a loop of magnetic field (where re is the radius of curvature of the field 
lines) can balance apressure gradient ~ P/rc, so the maximum pressure that can be contained 
in a magnetic fieldjis of order P ^ B^/87t. If the pressure is due to radiation, assuming an 
optical depth T > 1 so the diffusion approximation can be used, the maximum radiative energy 
flux is then of the order 

^r,max«cP,/T«C^. (13) 

In the range of validity of the assumptions made this has its maximum for an optical depth of 
order unity: /v.max ~ cB^ /%%. For a neutron star of radius R — 10^ cm and a field strength of 
lO'^ G, this gives Lr.max ~ lO'*^ erg/s, many orders of magnitude higher than the Eddington 
value Le- [This explains the enormous luminosities that can be reached in so-called magnetar 
outbursts, e.g. Hurley et al. 2005]. 

The Eddington argument considers only the radiative flux. Larger energy fluxes are possi- 
ble if energy is transported by other means, for example by convection. 

Since Le depends on opacity, it can happen that Le is lower in the atmosphere of a star than 
in its interior. A luminous star radiating near its (internal) Eddington rate will then blow off its 
atmosphere in a radiatively driven stellar wind; this happens for example in Wolf-Rayet stars. 
In the context of protostellar accretion, the opacity in the star forming cloud from which the 
protostar accretes is high due to atomic and molecular transitions. As a result, the radiation 
pressure from a massive (proto-)star, with a luminosity approaching ( |2.3| l, is able to clear the 
away the molecular cloud from which it formed. This is believed to set a limit on the mass 
that can be reached by a star formed in a molecular cloud. 



2,3.3 Neutron stars vs. black hole accreters 

In deriving the critical accretion rate, it was assumed that the gravitational energy liberated 
is emitted in the form of radiation. In the case of a black hole accreter, this is not necessary 
since mass can flow through the hole's horizon, taking with it all energy contained in it. 
Instead of being emitted as radiation, the energy adds to the mass of the hole. This becomes 
especially important at high accretion rates, M > Me (and in the ion supported accretion flows 



^ Depending on circumstances the actual maximum is less than this because a magnetically contained plasma 
tends to 'leak across' field lines through MHD instabilities. 
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discussed in section [T3]l. The parts of the flow close to the hole then become optically thick, 
the radiation stays trapped in the flow, and instead of producing luminosity gets swallowed 
by the hole. The accretion rate on a black hole can thus be arbitrarily large in principle (see 



also section 1 1 and chapter 10 in Frank et al, 2002). 

A neutron star cannot absorb this much energy (only a negligible amount is taken up by 
conduction of heat into its interior), so Mg is more relevant for neutron stars than for black 
holes. It is not clear to what extent it actually limits the possible accretion rate, however, 
since the limit was derived under the assumption of spherical symmetry. It is possible that 
accretion takes place in the form of an optically thick disk, while the energy released at the 
surface produces an outflow along the axis, increasing the maximum possible accretion rate 
(cf. discussion in[TT]). This has been proposed (e.g. King, 2004) as a possible conservative 
interpretation of the so-called ultraluminous X-ray sources (ULX), rare objects with lumi- 
nosities of 10^^ — 10^' erg/s. These are alternatively and more excitingly suggested to harbor 
intermediate mass black holes (above w 3QMq). 



3 Accretion with Angular Momentum 

When the accreting gas has a zonzero angular momentum with respect to the accreting object, 
it can not accrete directly. A new time scale then plays a role, the time scale for outward 
transport of angular momentum. Since this is in general much longer than the dynamical 
time scale, much of what was said about spherical accretion needs modification for accretion 
with angular momentum. 

Consider the accretion in a close binary consisting of a compact (white dwarf, neutron star 
or black hole) primary of mass Mi and a main sequence companion of mass M2. The mass 
ratio is defined as ^ = M2/M1 (note: in the literature q is just as often defined the other way 
around). 

If Ml and M2 orbit each other in a circular orbit and their separation is a, the orbital 
frequency Q is 

Q^ = G{Mi+M2)/a^. 

The accretion process is most easily described in a coordinate frame that corotates with this 
orbit, and with its origin in the center of mass. Matter that is stationary in this frame expe- 
riences an effective potential, the Roche potential (Ch. 4 in Frank, King and Raine, 2002), 
given by 

^ , , GM GM \. 2 

ri r2 2 

where ri 2 are the distances of point r to stars 1,2, and C7 the distance from the rotation axis 
(the axis through the center of mass, perpendicular to the orbit ). Matter that does not corotate 
experiences a very different force (due to the Coriolis force). The Roche potential is therefore 
useful only in a rather limited sense. For non-corotating gas intuition based on the Roche 
geometry can be misleading. Keeping in mind this Umitation, consider the equipotential sur- 



faces of ( 14 1. The surfaces of stars Mi .2, assumed to corotate with the orbit, are equipotential 



surfaces of ( 14 1. Near the centers of mass (at low values of 0r) they are unaffected by the 
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Other star, at higher they are distorted and at a critical value ^i the two parts of the surface 
touch. This is the critical Roche surface Si whose two parts are called the Roche lobes. 

Binaries lose angular momentum through gravitational radiation and a magnetic wind from 
the secondary (if it has a convective envelope). Through this loss the separation between the 
components decreases and both Roche lobes decrease in size. Mass transfer starts when M2 
fills its Roche lobe, and continues as long as the angular momentum loss from the system 
lasts. Mass transfer can also be due to expansion of the secondary in course of its evolution, 
and mass transfer can be a runaway process, depending on mass and internal structure of the 
secondary. This is a classical subject in the theory of binary stars, for an introduction see 
Warner (1995). 

A stream of gas then flows through the point of contact of the two parts of Si , the inner 
Lagrange point Li. If the force acting on it were derivable entirely from ( [T4| the gas would 
just fall in radially onto Mi . As soon as it moves however, it does not corotate any more and 
its orbit under the influence of the Coriolis force is different (Fig. [Til. 




Fig. 1 Roche geometry for q = 0.2, with free particle orbit from Li (as seen in a frame corotating with the 
orbit). Dashed: circularization radius. 



Since the gas at Li is usually very cold compared with the virial temperature, the veloc- 
ity it acquires exceeds the sound speed already after moving a small distance from Li. The 
flow into the Roche lobe of Mi is therefore highly supersonic. Such hypersonic flow is essen- 
tially ballistic, that is, the stream flows approximately along the path taken by freely falling 
particles. 

Though the gas stream on the whole follows an path close to that of a free particle, a strong 
shock develops at the point where the path intersects itselnl After this, the gas settles into a 
ring, into which the stream continues to feed mass. If the mass ratio q is not too small this ring 



In practice shocks already develop shortly after passing the pericenter at M] , when the gas is decelerated 
again. Supersonic flows that are decelerated, by whatever means, in general develop shocks (e.g. Courant and 
Friedrichs 1948, Massey, 1968). The effect can be seen in action in the movie published in Rozyczka and 
Spruit, 1993. 
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forms fairly close to Mi . An approximate value for its radius is found by noting that near Mi 
the tidal force due to the secondary is small, so that the angular momentum of the gas with 
respect to Mi is approximately conserved. If the gas continues to conserve angular momentum 
while dissipating energy, it settles into the minimum energy orbit with the specific angular 
momentum j of the incoming stream. The value of j is found by a simple integration of the 
orbit starting at Li and measuring j at some point near pericenter. The radius of the orbit, 
the nominal circularization radius r^ is then defined through (GMirc)' " = j. In units of the 
orbital separation a, r^ and the distance tli from Mi to Li are functions of the mass ratio only. 
As an example for q = 0.2, rn ~ 0.66fl and the circularization radius r^ w 0.16fl. In practice 
the ring forms somewhat outside r^ because there is some angular momentum redistribution 
in the shocks that form at the impact of the stream on the ring. The evolution of the ring 
depends critically on nature and strength of the angular momentum transport processes. If 
sufficient 'viscosity' is present it spreads inward and outward to form a disk. 

At the point of impact of the stream on the disk the energy dissipated is a significant 
fraction of the orbital kinetic energy, hence the gas heats up to a significant fraction of the 
virial temperature. For a typical system with Mi ~ IMq,, M2 = Q.2Mq and an orbital period 
of 2 hrs, the observed size of the disk (e.g. Wood et al. 1989b, Rutten et al. 1992) r^/a « 0.3, 
the orbital velocity at r^ about 900 km/s, and the virial temperature at r^ is ss 10^ K. The 
actual temperatures at the impact point are much lower, due to rapid cooling of the shocked 
gas. Nevertheless the impact gives rise to a prominent 'hot spot' in many systems, and an 
overall heating of the outermost part of the disk. 



4 Thin disks: properties 

4.1 Flow in a cool disk is supersonic 

Ignoring viscosity, the equation of motion in the potential of a point mass is 

dy „ 1„ GM^ 

+ vVv = VP-— ^f, (15) 

dt p r^ 

where f is a unit vector in the spherical radial direction r. To compare the order of magnitude 
of the terms, choose a position ro in the disk, and choose as typical time and velocity scales 
the orbital time scale H^ — (jqIGMY''^ and velocity i2oro. For simplicity of teh argument, 
assume a fixed temperature T . The pressure gradient term is then 

1 ^ 

-VP== — rvinP. 

p M 

In terms of the dimensionless quantities 

r^r/m, i)=: v/(i2oro), 
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the equation of motion becomes 

^+v-Vv= VlnP-^f. (16) 

at Tvir F 

All terms and quantities in this equation are of order unity by the assumptions made, except 
the pressure gradient term which has the coefficient T/Tyir in front. If cooling is important, so 
that r/Tvir <C 1, the pressure term is negligible to first approximation. Equivalent statements 
are also that the gas moves hypersonically on nearly Keplerian orbits, and that the disk is thin, 
as is shown next. 



4.2 Disk thickness 

The thickness of the disk is found by considering its equilibrium in the direction perpendicu- 
lar to the disk plane. In an axisymmetric disk, using cylindrical coordinates (CJ, ^ , z), consider 
the forces at a point ro (OJ, 0,0) in the midplane, in a frame rotating with the Kepler rate Qq 
at that point. The gravitational acceleration —GM/r^f balances the centrifugal acceleration 
Qq us at this point, but not at some distance z above it because gravity and centrifugal acceler- 
ation work in different directions. Expanding both accelerations near ro, one finds a residual 
acceleration toward the midplane of magnitude 

gz = -^Iz- 

Assuming an isothermal gas at temperature T, the condition for equilibrium in the z- direction 
under this acceleration yields a hydrostatic density distribution 

p=Po(OT)expf-^j. 

H{U5), called the scale height of the disk or simply the disk thickness, is given in terms of the 
isothermal sound speed c\ — {MT /jxyi'^ by 

H = cJQq. (17) 

Define 5 = H /r, the aspect ratio of the disk; it can be expressed in several equivalent ways: 

«.«.^.>.m"; OS) 

r Llr Ma \T^ir J 
where Ma is the Mach number of the orbital motion. 
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4.3 Viscous spreading 

The shear flow between neighboring Kepler orbits in the disk causes friction if some form 
of viscosity is present. The frictional torque is equivalent to exchange of angular momen- 
tum between these orbits. But since the orbits are close to Keplerian, a change in angular 
momentum of a ring of gas also means it must change its distance from the central mass. If 
the angular momentum is increased, the ring moves to a larger radius. In a thin disk angular 
momentum transport (more precisely a nonzero divergence of the angular momentum flux) 
therefore automatically implies redistribution of mass in the disk. 

A simple example (Liist 1952, see also Lynden-Bell and Pringle 1974) is a narrow ring of 
gas at some distance tq. If at f = this ring is released to evolve under the viscous torques, 
one finds that it first spreads into an asymmetric hump. The hump quickly spreads inward 
onto the central object while a long tail spreads slowly outwards to large distances. As f ^ oo 
almost all the mass of the ring accretes onto the center, while a vanishingly small fraction of 
the gas carries almost all the angular momentum to infinity. 

As a result of this asymmetric behavior essentially all the mass of a disk can accrete even 
if the total angular momentum of the disk is conserved. In practice, however, there is often an 
external torque removing angular momentum from the disk: when the disk results from mass 
transfer in a binary system. The tidal forces exerted by the companion star take up angular 
momentum from the outer parts of the disk, limiting its outward spreading. 



4.4 Observational evidence of disk viscosity 

Evidence for the strength of the angular momentum transport processes in disks comes 
from observations of variability time scales. This evidence is not good enough to determine 
whether the processes really behave in the same way as viscosity, but if this is assumed, 
estimates can be made of its magnitude. 

Observations of Cataclysmic Variables (CV) give the most detailed information. These 
are binaries with white dwarf (WD) primaries and (usually) main sequence companions (for 
reviews see Meyer-Hofmeister and Ritter 1993, 1995, Warner 1995). A subclass of these sys- 
tems, the Dwarf Novae, show semiregular outbursts. In the currently most developed theory, 
these outbursts are due to an instability in the disk (Smak 1971, Meyer and Meyer-Hofmeister 
1981, King 1995, Hameury et al. 1998). The outbursts are episodes of enhanced accretion of 
mass from the disk onto the primary, involving a significant part of the whole disk. The de- 
cay time of the burst is thus a measure of the viscous time scale of the disk (the quantitative 
details depend on the model, see Cannizzo et al. 1988, Hameury et al. 1998): 

fvisc — ''d/^i 

where r^ is the size of the disk, ^ lO'" cm for a CV. With decay times on the order of 
days, this yields viscosities of the order lO'^ cm^/s, some 14 orders of magnitude above the 
microscopic ('molecular') viscosity of the gas. 
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Other evidence comes from the inferred time scale on which disks around protostars dis- 
appear, which is of the order of 10^ years (e.g. Strom et al., 1993). 



4.5 a-parameterization 

Several processes have been proposed to account for these short time scales and the large ap- 
parent viscosities inferred from them. One of these is that accretion does not in fact take place 
through a viscous-like process as described above at all, but results from angular momentum 
loss through a magnetic wind driven from the disk surface (Bisnovatyi-Kogan & Ruzmaikin 
1976, Blandford 1976), much in the way sun-like stars spin down by angular momentum loss 
through their stellar winds. The extent to which this plays a role in accretion disks is still 
uncertain. It would be a 'quiet' form of accretion, since it can do without energy dissipating 
processes like viscosity. It has been proposed as the explanation for the low ratio of X-ray 
luminosity to jet power in many radio sources (see Migliari and Fender 2006 and references 
therein). This low ratio, however, is also plausibly attributed to the low efficiency with which 



X-rays are produced in the 'ion supported accretion flows' (see section 13 below) which are 
expected to be the source of the jet outflows observed from X-ray binaries and AGN. 

The most quantitatively developed mechanism for angular momentum transport is a form 
of magnetic viscosity (anticipated already in Shakura and Sunyaev 1973). This is discussed 
below in section [8T| It requires the accreting plasma to be sufficiently electrically conducting. 
This is often the case, but not always: it is questionable for example in the cool outer parts 
of protostellar disks. Other mechanisms thus still play a role in the discussion, for example 
spiral shocks (Spruit et al. 1987) and self-gravitating instabilities (Paczyriski 1978, Gammie 
1997). 

In order to compare the viscosities of disks in objects with (widely) different sizes and 
physical conditions, introduce a dimensionless viscosity a: 

v = a|, (19) 

where c; is the isothermal sound speed as before. The quantity a was introduced by Shakura 
and Sunyaev (1973) as a way of parametrizing our ignorance of the angular momentum trans- 
port process, in a way that allows comparison between systems of very different size and 
physical origin (Their definition of a differs a bit, by a constant factor of order unity). 



4.6 Causality limit on turbulent viscosity 

How large can the value of a be, on theoretical grounds? As a simple model, let's assume 
that the shear flow between Kepler orbits is unstable to the same kind of shear instabilities 
found for flows in tubes, channels, near walls and in jets. These instabilities occur so ubiq- 
uitously that the fluid mechanics community considers them an automatic consequence of a 
high Reynolds number: 
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LV 
Re= — 
V 

where L and V are characteristic length and velocity scales of the flow. If this number exceeds 
about 1000 (for some forms of instability much less), instability and turbulence are generally 
observed. It has been argued (e.g. Zel'dovich 1981) that for this reason hydrodynamic turbu- 
lence is the cause of disk viscosity. Let's look at the consequences of this assumption. If an 
eddy of radial length scale I develops due to shear instability, it will rotate at a rate given by 
the rate of shear in the flow, C7, in a Keplerian disk 

or 2 

The velocity amplitude of the eddy is V = al, and a field of such eddies would produce a 
turbulent viscosity of the order (leaving out numerical factors of order unity) 

V,urb«^^i2. (20) 

In compressible flows, there is a maximum to the size of the eddy set by causality consider- 
ations. The force that allows an instability to form a coherently overturning eddy is the pres- 
sure, which transports information about the flow at the speed of sound. The eddies formed 
by a shear instability can therefore not rotate faster than Ci, hence their size does not exceed 



c\/o « H (eq. 18 1. At the same time, the largest eddies formed also have the largest con- 
tribution to the exchange of angular momentum. Thus we should expect that the turbulent 
viscosity is given by eddies with size of the order H: 

v<H^-n, 

or 

a < 1. 

The dimensionless number a can thus be interpreted as the effective viscosity in units of the 
maximum value expected if a disk were hydrodynamically turbulent. 



4,6.1 Large scale vortices? 

The small size of hydrodynamic eddies expected, L < H ^ c^/Q, is due to the high Mach 
number of the flow in an accretion disk; this makes a disk behave like a very compress- 
ible fluid. Attempts to construct large scale vortices L^ H in disks using incompressible 
fluid analogies continue to be made, both analytically and experimentally. Expansions in disk 
thickness as a small parameter then suggest themselves, in analogy with large scale flows 
such as weather systems and tropical storms in the Earth's atmosphere. The size of these 
systems is large compared with the height of the Earth's atmosphere, and expansions making 
use of this can be effective. This tempting analogy is misleading for accretion disks, however, 
since in contrast with the atmosphere all flows with a horizontal scale exceeding the vertical 
thickness are supersonic, and the use of incompressible fluid models meaningless. 
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4.7 Hydrodynamic turbulence? 

Does hydrodynamical turbulence along these lines actually exist in disks? In the astrophysical 
community, a consensus has developed that it does not, certainly not along the simple lines 
suggested by a Reynolds' number argument. An intuitive argument is that the flow in a cool 
disk is close to Kepler orbits, and these are quite stable. Numerical simulations of disk-like 
flows with different rotation profiles do not show the expected shear flow instabilities in cases 
where the angular momentum increases outward, while turbulence put in by hand as initial 
condition decays in such simulations (Hawley et al. 1999). 

Analytical work on the problem has not been able to demonstrate instability in this case 
either (for recent work and references see Lesur & Longaretti 2005, Rincon et al. 2007). 
Existing proposals for disk turbulence involve ad-hoc assumptions about the existence of 
hydrodynamic instabilities (e.g. Dubrulle 1992). However, a subtle form of hydrodynamic 
angular momentum transport has been identified more recently (Lesur & Papaloizou 2010); it 
depends on vertical stratification of the disk as well as the presence of a convectively unstable 
radial gradient. 

A laboratory analogy is the rotating Couette flow, an experiment with water between dif- 
ferentially rotating cylinders. Recent such experiments have demonstrated that turbulence is 
absent in cases where angular momentum increases outward (as in Keplerian rotation) at 
Reynolds' numbers as high as 10^ (the Princeton Couette experiment, Ji et al. 2006). 

In view of these negative results, the popular mechanism for angular momentum transport 



in accretion disk has become magnetic: 'MRI turbulence' (section 8.1 1. But beware: in the 
fluid mechanics community it is considered crackpot to suggest that at Reynolds numbers like 
our 10^^ a flow (no matter which or where) can be anything but 'fiercely turbulent'. If your 
career depends on being friends with this community, it would be wise to avoid discussions 
about accretion disks or Couette experiments. 



5 Thin Disks: equations 



Consider a thin (= cool, nearly Keplerian, cf. section 4.2 1 disk, axisymmetric but not station 



ary. Using cylindrical coordinates {r,(p,z), (note that we have changed notation from CT to r 



compared with |4.2[ l we define the surface density E of the disk as 

pdz^lHopo, (21) 



where po, Hq are the density and scaleheight at the midplane. The approximate sign is used 
to indicate that the coefficient in front of H in the last expression actually depends on details 
of the vertical structure of the disk. Conservation of mass, in terms of L is described by 

^(r£) + ^(ri:vv)=0. (22) 
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(derived by integrating the continuity equation over z). Since the disk is axisymmetric and 
nearly Keplerian, the radial equation of motion reduces to 

vl - GM/r. (23) 

The ^-equation of motion is 

dt or r 

where F^ is the azimuthal component of the viscous force. By integrating this over height z 
and using ( [22] i, one gets an equation for the angular momentum balance: 

^{rLar^) + ^{rLv,ar^) = ^{Sr'^), (25) 

dt dr dr dr 

where Q = v^/r, and 

S= [ pvdz^Lv. (26) 



The second approximate equality in ([26]) holds if v can be considered independent of z- 



The right hand side of ( 25 i is the divergence of the viscous angular momentum flux, and is 
derived most easily with a physical argument, as described in, e.g. Pringle (1981) or Frank et 
al. (2002f] 

Assume now that v can be taken constant with height. For an isothermal disk {T indepen- 
dent of z), this is equivalent to taking the viscosity parameter a independent of z- As long 
as we are not sure what causes the viscosity this is a reasonable simplification. Note, how- 
ever, that recent numerical simulations of magnetic turbulence suggest that the effective a, 
and the rate of viscous dissipation per unit mass, are higher near the disk surface than near 



the midplane (c.f. section^. While eq. (25 1 is still valid for rotation rates Q deviating from 
Keplerian (only the integration over disk thickness must be justifiable), we now use the fact 
that i2 ~ r^^". Then eqs. (22|25 1 can be combined into a single equation for Z: 



.¥=3^[r'/4(vr.'/2)]. (27) 



dt dr d 



r 



Under the same assumptions, eq. ( 24 1 yields the mass flux M at any point in the disk: 



' If you prefer a more formal derivation, the fastest way is to consult Landau and Lifshitz (1959) chapter 15 
(hereafter LL). Noting that the divergence of the flow vanishes for a thin axisymmetric disk, the viscous stress 
a becomes (LL eq. 15.3) 

(dvj dvi, 
^r + 1^ 
dxk dxi 

where r\ = pv. This can be written in cylindrical or spherical coordinates using LL eqs. (15.15-15.18). The 

viscous force is 

(9c7,T 1 dr] „-, 

dxk T] dxk 

Writing the Laplacian in cylindrical coordinates, the viscous torque is then computed from the (^-component 

of the viscous force by multiplying by r, and is then integrated over z- 
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M^-2nrEvr = 6nr^/^ — (vEr^^^). (28) 

dr 



Eq. ( 27 1 is the standard form of the thin disk diffusion equation. An important conclusion 
from this equation is: in the thin disk Hmit, all the physics which influences the time depen- 
dent behavior of the disk enters through one quantity only; the viscosity v. This is the main 
attraction of the thin disk approximation. 



5.1 Steady thin disks 

In a steady disk (d/dt = 0) the mass flux M is constant through the disk and equal to the 
accretion rate onto the central object. From ([28|) we get the surface density distribution: 



vZ 



3;r 


1- 


AT] 



(29) 



where r, is the inner radius of the disk and )3 is a parameter appearing through the integration 
constant. It is related to the flux of angular momentum Fj through the disk: 

Fj = -MpQtrf, (30) 

where X2, is the Kepler rotation rate at the inner edge of the disk. If the disk accretes onto an 
object with a rotation rate i2* less than i2,, the rotation rate Q (r) as a function of distance r 
jumps from X2* close to the star to the Kepler rate QK{r) in the disk. It then has a maximum 
at some distance close to the star (distance of the order of the disk thickness H). At this 
point the viscous stress (proportional to drQ) vanishes and the angular momentum flux is 



just the amount carried by the accretion flow. In (30 1 this corresponds to j3 = 1, independent 
of 12* (Shakura and Sunyaev, 1973, Lynden-Bell and Pringle, 1974). This is referred to as the 
standard or 'accreting case' . The angular momentum flux (equal to the torque on the accreting 
star), is then inward, causing the rotation rate of the star to increase (spinup). 

For stars rotating near their maximum rate (X2* « X2,) and for accretion onto magneto- 
spheres, which can rotate faster than the disk, the situation is different (Sunyaev and Shakura 
1977, Popham and Narayan 1991, Paczyhski 1991, Bisnovatyi-Kogan 1993, Rappaport et al. 
2004). If the inner edge of the disk is at the corotation radius r^o, defined by 

X2K(rco)=i2*, (31) 

the viscous stress cannot be assumed to vanish there, and the thin disk approximation does not 
give a unique answer for the angular momentum flux parameter j3 . Its value is then determined 
by details of the hydrodynamics at r^o, which need to be investigated separately. Depending 
on the outcome of this investigation, values varying from 1 (spinup, standard accreting case) 
to negative values (spindown) are possible. As ([30]l shows, the surface density at the inner 
edge of the disk depends sensitively on the value of /3 (see also Fig. l2|i. This plays a role in 
the cyclic accretion process discussed in the next subsection. 
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Fig. 2 Surface density vZ of a tliin disc as a function of distance from the corotation radius r^ , for a steady, 
thin viscous disc. Left: steady accretion at a fixed accretion rate m, for inner edge of the disc at corotation. /3 
measures the angular momentum flux, j8 = — 1 corresponding to the standard case of accretion on to a slowly 
rotation object. For ji < the angular momentum flux is outward (spindown of the star). Right: 'quiescent 
disc' solutions with m = and a steady outward angular momentum flux due to a torque / applied at the inner 
edge. 



5.2 Magnetospheres, 'pwpellering' and 'dead disks' 



Stars with magnetospheres, instead of spinning up by accretion of angular momentum from 
the disk, can actually spin down by interaction with the disk, even while accretion is still 
going on. The thin disk approximation above covers this case as well (cf. figure l2]i. The 
surface density distribution is then of the form (29 1, but with j3 < (see also Spruit and Taam 
1993, Rappaport et al. 2004). The angular momentum flux through the disk is outward, and 
the accreting star spins down. This is possible even when the interaction between the disk 
and the magnetosphere takes place only at the inner edge of the disk. Magnetic torques due 
interaction between the disk and the magnetosphere may exist at larger distances in the disk 
as well, but are not necessary for creating an outward angular momentum flux. Numerical 
simulations of disk-magnetosphere interaction (Miller and Stone 1997, for recent work see 
Long et al. 2008 and references therein) give an interesting view of how such interaction may 
take place, presenting a picture that is very differently from what is assumed in the previous 
'standard' models. Among other things, they show the interaction region to be quite narrow. 
Rapidly rotating magnetospheres play a role in some CVs and X-ray binaries, and probably 
also in protostellar accretion disks. When the rotation velocity of the magnetosphere is larger 
than the Kepler velocity at the inner edge of the disk (i.e. when r; > rco), the literature often 
assumes that the mass transfered from the secondary must be 'flung out' of the binary system. 
This idea is idea called 'propellering' (lUarionov and Sunyaev 1975). The term is then used 
as synonymous with the condition ri > r^o. While a process like this is likely to happen when 
the rotation rate of the star is sufficiently large (the CV AE Aqr being an example) it is not 
necessary. It is not possible either, when the difference in rotation velocity is too small to put 
the accreting mass on an escape orbit. This has been realized early on (Sunyaev and Shakura 
1977). Instead of being flung out, accretion is halted and mass accumulates in the disk. The 
buildup of mass in the disk then leads to a following episode of accretion. In this phase the 
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Fig. 3 Cyclic accretion due to interaction between a magnetosphere and a disk with inner edge near coro- 
tation. Solid: inner edge radius of the disk in units of the corotation radius; dashed: accretion rate (from 
D'Angelo and Spruit 2010). 



magnetosphere is compressed such that n < rco- Instead of propellering, the result is cyclic 
accretion. For quantitative models see D'Angelo and Spruit (2010), an example is shown in 
FigurelS] The phase with r^ > rco, during which accretion is halted is called 'dead disk' phase 
in Sunyaev and Shakura (1977). 



5.3 Disk Temperature 



In this section I assume accretion onto not-too-rapidly rotating objects, so that j3 = 1 (eq. 29 1 



The surface temperature of the disk, which determines how much energy it loses by radiation, 
is governed by the energy dissipation rate in the disk, which in turn is given by the accretion 
rate. From the first law of thermodynamics we have 



dS 
df 



(32) 



where S is the entropy per unit mass, F the heat flux (including radiation and convection), 
and Qy the viscous dissipation rate. For thin disks, where the advection of internal energy 
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(or entropy) can be neglected, and for changes which happen on time scales longer than the 
dynamical time Q^^, the left hand side is small compared with the terms on the right hand 
side. Integrating over z, the divergence turns into a surface term and we get 

2cJrr/ = /_"evdz, (33) 

where Tg is the surface temperature of the disk, Oi- is Stefan-Boltzmann's radiation constant 
<7r — a^c/A, and the factor 2 comes about because the disk has 2 radiating surfaces (assumed 
to radiate like black bodies). Thus the energy balance is local (for such slow changes): what 
is generated by viscous dissipation inside the disk at any radius r is also radiated away from 
the surface at that position. The viscous dissipation rate is equal to Q^, = aijdvi/dxj, where 
Oij is the viscous stress tensor (see footnote in section 5), and this works ourlto be 

Q,=9/AQ^vp. (34) 



Eq. ( 33 I, using ( 29 1 then gives the surface temperature in terms of the accretion rate: 

'r,-\l/2l 



4 9^7 ^ GM3M 



1 

r 



(35) 



This shows that the surface temperature of the disk, at a given distance r from a steady ac- 
creter, depends only on the product MM, and not on the highly uncertain value of the viscosity. 
For r :^ rj we have 

T, ~ r-^/\ (36) 

These considerations only tells us about the surface temperature. The internal temperature 
in the disk is quite different, and depends on the mechanism transporting energy to the sur- 
face. Because it is the internal temperature that determines the disk thickness H (and probably 
also the viscosity), this transport needs to be considered in some detail for realistic disk mod- 
els. This involves a calculation of the vertical structure of the disk. Because of the local (in 
r) nature of the balance between dissipation and energy loss, such calculations can be done 
as a grid of models in r, without having to deal with exchange of energy between neighbor- 
ing models. Schemes borrowed from stellar structure computations are used (e.g. Meyer and 
Meyer-Hofmeister 1982, Cannizzo et al. 1988). 

An approximation to the temperature in the disk can be found when a number of additional 
assumptions is made. As in stellar interiors, the energy transport is radiative rather than con- 
vective at high temperatures. Assuming local thermodynamic equilibrium (LTE, e.g. Rybicki 
and Lightman 1979), the temperature structure of a radiative atmosphere is given, in the so- 
called Eddington approximation of radiative transfer (not to be confused with the Eddington 
limit) by: 

l.,r'.5K ,37) 

The boundary condition that there is no incident flux from outside the atmosphere yields the 
approximate condition 



using, e.g. LL eq. 16.3 
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aj\x^2/3)^F, 
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(38) 



where T = /_°° Kpdz is the optical depth at geometrical depth z, and F the energy flux through 
the atmosphere. Assuming that most of the heat is generated near the midplane (which is the 
case if v is constant with height), F is approximately constant with height and equal to O-iT^, 



given by ( 35 i. Eq ( 37 1 then yields 



aj'^-{T+-)F. 



(39) 



Approximating the opacity K as constant with z, the optical depth at the midplane is T 
kT,/2. If t ^ 1, the temperature at the midplane is then: 



64 



(40) 



With the equation of state (J2|l, valid when radiation pressure is small, we find for the disk 
thickness, using (29 1: 



- = (^/Ai)2/5((j,64/37r2)"'/'''((f/a)i/iO(GM)-^/20rV20(/M)i/5 

= 5lQ-^a-'l'''rf\M/MQ)-y^\fM,^Yl\ [P,. ^ P) (41) 

where re = r/{lQ^ cm), Mig = M/(10''^g/s), and 

/=l-(r,/r)^/2. 

From this we conclude: i) that the disk is thin in X-ray binaries, H/r < 0.01, ii) the disk 
thickness is relatively insensitive to the parameters, especially a, K and r. It must be stressed, 
however, that this depends fairly strongly on the assumption that the energy is dissipated in 
the disk interior. If the dissipation takes place close to the surface, such as in some magnetic 
reconnection models, the internal disk temperature will be much closer to the surface tem- 
perature (Haardt et al. 1994, Di Matteo et al. 1999a and references therein). The midplane 



temperature and H are even smaller in such disks than calculated from (41 1 



5.4 A subtlety: viscous vs. gravitational energy release 



The viscous dissipation rate per unit area of the disk, Wy = {9/A)Q^vE [cf. eq. 35 1] can be 
compared with the local rate Wq at which gravitational energy is liberated in the accretion 
flow. Since half the gravitational energy stays in the flow as orbital motion, we have 



Wg 



1 GMM 

27tr 2r^ 



(42) 



so that, in a steady thin disk accreting on a slowly rotating object (using eq. 29 1: 
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H'v/Wg - 3/ = 3 [1 - (n/r) ^/2] . (43) 

At large distances from the inner edge, the dissipation rate is 3 times larger than the rate 
of gravitational energy release. This may seem odd, but becomes understandable when it is 
realized that there is a significant flux of energy through the disk associated with the viscous 
stress^] Integrating the viscous energy dissipation over the whole disk, one finds 



/: 



GMM 

2nrW^Ar=^- , (44) 

2ri 



as expected. That is, globally but not locally, half of the gravitational energy is radiated from 
the disk while the other half remains in the orbital kinetic energy of the material accreting on 
the star. 

What happens to this remaining orbital energy depends on the nature of the accreting 
object. If the object is a slowly rotating black hole, the orbital energy is just swallowed by the 
hole. If it has a solid surface, the orbiting gas slows down until it corotates with the surface, 
dissipating the orbital energy into heat in a boundary layer. Unless the surface rotates close 
to the orbital rate ('breakup'), the energy released in this way is of the same order as the 
total energy released in the accretion disk. The properties of this boundary layer are therefore 



critical for accretion onto neutron stars and white dwarfs. See also section 9.1 and Inogamov 
and Sunyaev (1999). 



5.5 Radiation pressure dominated disks 

In the inner regions of disks in XRB, the radiation pressure can dominate over the gas pres- 
sure, which results in a different expression for the disk thickness. The total pressure P is 



P = P,.+Pf,= -ar+Pg. (45) 



Defining a 'total sound speed' by cl = P/p the relation between temperature and disk thick- 
ness (flTli stin holds, so Ct = Q.H. For Pr > Pg we get from (pOl), with (|35]) and T > 1: 



(where the rather approximate relation Z = 2//po has been used). Thus, 

"i^^^fM^'-ft, (46) 

R StzcR-' 2-' Me 

where R is the stellar radius and Me the Eddington rate for this radius. It follows that the disk 
becomes thick near the star, if the accretion rate is near Eddington (though this is mitigated 



^ See Landau & Lifshitz section 16 
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somewhat by the factor /). Accretion near the Eddington Hmit is evidently not geometrically 
thin any more. In addition, other processes such as angular momentum loss from the disk by 
'photon drag' in the dense radiation field of the accreting star have to be taken into account. 



6 Time scales in thin disks 

A number of time scales play a role in the behavior of disks. In a thin disk they differ by 
powers of the large factor r/H. The shortest of these is the dynamical time scale t^: 

fd = ^K"' = {GM/r-y^l^. (47) 

The time scale for radial drift through the disk over a distance of order r is the viscous time 
scale: 



(using (28 and (29 1, valid for steady accretion). This is of the order [rjH)^ longer than the 
dynamical time scale, except near the inner edge of the disk, where / 4 0, andF this time scale 
formally drops to zero as a consequence of the thin disk approximations. If more physics were 
included, the viscous time scale would stay above the dynamical time scale. 

Intermediate are thermal time scales. If fit is the thermal energy content (enthalpy) of the 
disk per unit of surface area, and Wy = (9/4)i2^vZ the heating rate by viscous dissipation, 
we can define a heating time scale: 

fh=£t/Wv. (49) 

In the same way, a cooling time scale is defined by the energy content and the radiative loss 
rate: 

rc=£t/(2cj,7;4). (50) 

For a thin disk, the two are equal since the viscous energy dissipation is locally balanced by 
radiation from the two disk surfaces. [In thick disks (ADAFs), this balance does not hold, 
since the advection of heat with the accretion flow is not negligible. In this case fc > fh (see 
sectionfTS])]. Thus, we can replace both time scales by a single thermal time scale ft, and find, 
wifli(|34|: 

1 r jp 

h = —\ ^dz, (51) 

where the enthalpy of an ideal gas of constant ratio of specific heats 7 has been used. Leaving 
out numerical factors of order unity, this yields 

,.^. (52) 

That is, the thermal time scale of the disk is independent of most of the disk properties and 
of the order \ja times longer than the dynamical time scale. This independence is a conse- 
quence of the a-parameterization used. If a is not a constant, but dependent on disk temper- 
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ature for example, the dependence of the thermal time scale on disk properties will become 
apparent again. 

If, as seems likely from observations, a is generally < 1, we have in thin disks the ordering 
of time scales; 

fv > ft > fd- (53) 



7 Comparison with CV observations 

The number of meaningful quantitative tests between the theory of disks and observations 
is somewhat limited since in the absence of a theory for v, it is a bit meagre on predictive 
power The most detailed information perhaps comes from modeling of CV outbursts. 

Two simple tests are possible (nearly) independently of v. These are the prediction that 



the disk is geometrically quite thin (eq. 41 1 and the prediction that the surface temperature 
Ts ~ r"^'"* in a steady disk. The latter can be tested in a subclass of the CV's that do not show 
outbursts, the nova-like systems, which are believed to be approximately steady accreters. 
If such a system is also eclipsing, eclipse mapping techniques can be used to derive the 
brightness distribution with r in the disk (Home, 1985, 1993). If this is done in a number of 
colors so that bolometric corrections can be made, the results (e.g. Rutten et al. 1992) show 
in general a fair agreement with the r^^'"* prediction. Two deviations occur: i) a few systems 
show significantly flatter distributions than predicted, and ii) most systems show a 'hump' 
near the outer edge of the disk. The latter deviation is easily explained, since we have not 
taken into account that the impact of the stream additionally heats the outer edge of the disk. 
Though not important for the total light from the disk, it is an important local contribution 
near the edge. 

Eclipse mapping of Dwarf Novae in quiescence gives a quite different picture. Here, the 
inferred surface temperature profile is often nearly flat (e.g. Wood et al. 1989a, 1992). This 
is understandable however since in quiescence the mass flux depends strongly on r. In the 
inner parts of the disk it is small, near the outer edge it is close to its average value. With 



eq. (35 1, this yields a flatter Ts{r). The lack of light from the inner disk is compensated 
during the outburst, when the accretion rate in the inner disk is higher than average (see 
Mineshige and Wood 1989 for a more detailed comparison). The effect is also seen in the 
2-dimensional hydrodynamic simulations of accretion in a binary by Rozyczka and Spruit 
(1993). These simulations show an outburst during which the accretion in the inner disk is 
enhanced, between two episodes in which mass accumulates in the outer disk. 



7.1 Comparison with LMXB observations: irradiated disks 

In low mass X-ray binaries a complication arises because of the much higher luminosity of 
the accreting object. Since a neutron star is roughly 1000 times smaller than a white dwarf, it 
produces 1000 times more luminosity for a given accretion rate. 
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Irradiation of the disk by the central source (the accreting star plus inner disk) leads to a 



different surface temperature than predicted by (35 i. The central source radiates nearly the 
total accretion luminosity GMM/R (assuming sub-Eddington accretion, see section 2). If the 
disk is concave, it will intercept some of this luminosity. If the central source is approximated 
as a point source the irradiating flux on the disk surface is 

F -'-^^ (54) 

where e is the angle between the disk surface and the direction from a point on the disk 
surface to the central source: 

e^dH/dr-H/r. (55) 

The disk is concave if £ is positive. We have 



F 3fR' 



where F is the flux generated internally in the disk, given by ( 35 1. On average, the angle e is 
of the order of the aspect ratio 5 = H/r. With / ~ 1, and our fiducial value 5 sa 5 10^^, we 
find that irradiation in LMXB dominates for r ^ 10^ cm. This is compatible with observations 
(for reviews see van Paradijs and McClintock 1993), which show that the optical and UV are 
dominated by reprocessed radiation from the innermost regions. 

When irradiation by an external source is included in the thin disk model, the surface 



boundary condition of the radiative transfer problem, equation ( 38 1 becomes 

(j,7;4=F + (l-fl)F„, (56) 

where a is the X-ray albedo of the surface, that is, 1 — a is the fraction of the incident flux that 
is absorbed in the optically thick layers of the disk (photons absorbed higher up only serve to 
heat up the corona of the diskrl The surface temperature T^ increases in order to compensate 
for the additional incident heat flux. The magnitude of the incident flux is sensitive to the 
assumed disk shape H{r), as well as on the assumed shape (plane or spherical, for example) 
of the central X-ray emitting region. 

The disk thickness depends on temperature, and thereby also on the irradiation. It turns out, 
however, that this dependence on the irradiating flux is small, if the disk is optically thick and 



the energy transport is by radiation (Lyutyi and Sunyaev 1976). To see this, integrate (37i 



with the modified boundary condition (56 1. This yields 



(7,r'* = ^F(T+^) + (l-fl)F„. (57) 

Thus the irradiation adds an additive constant to T'^{z). At the midplane, this constant has 
much less effect than at the surface. For the midplane temperature and the disk thickness to 
be affected significantly, it is necessary that 



* Incorrect derivations exist in tlie literature, in wliich the effect of irradiation is treated like an energy flux 
added to the internal viscous dissipation inside the disk rather than incident on the surface. 
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Fin/F ;> T. (58) 

The reason for this weak dependence of the midplane conditions on irradiation is the same 
as in radiative envelopes of stars, which are also insensitive to the surface boundary condi- 
tion. The situation is different for convective disks. As in fully convective stars, the adiabatic 
stratification then causes the conditions at the midplane to depend much more directly on 
the surface temperature. The outer parts of the disks in LMXB with wide orbits are in fact 
convective, hence their thickness more directly affected by irradiation. This leads to the pos- 
sibility of shadowing, with irradiated, vertically extended regions blocking irradiation of the 
disk outside. This plays an observable role in protostellar disks (Dullemond et al. 2001). 



7.2 Observational evidence of disk thickness 

From the paucity of sources in which eclipses of the central source by the companion are 
observed one deduces that the companion is barely or not at all visible as seen from the 
inner disk. Apparently some parts of the disk are much thicker than expected from the above 
arguments. This is consistent with the observation that a characteristic modulation of the 
optical light curve indicative of irradiation of the secondary's surface by the X-rays is not 
very strong in LMXB (an exception being Her X-1, which has an atypically large companion). 
The place of the eclipsing systems is taken by the so-called 'Accretion Disk Corona' (ADC) 
systems, where shallow eclipses of a rather extended X-ray source are seen instead of the 
expected sharp eclipses of the inner disk (for reviews of the observations, see Lewin et al. 
1995). The conclusion is that there is an extended X-ray scattering 'corona' above the disk. 
It scatters a few per cent of the X-ray luminosity. 

What causes this corona and the large inferred thickness of the disk? The thickness ex- 
pected from disk theory is a rather stable small number. To 'suspend' matter at the inferred 
height of the disk forces are needed that are much larger than the pressure forces available 
in an optically thick disk. A thermally driven wind, produced by X-ray heating of the disk 
surface, has been invoked (Begelman et al. 1983, Schandl and Meyer 1994, 1997). For other 
explanations, see van Paradijs and McClintock (1995). Perhaps a magnetically driven wind 
from the disk (e.g. Bisnovatyi-Kogan & Ruzmaikin 1976, Blandford & Payne 1982), such 
as inferred for protostellar objects (cf. Lee et al. 2000) can explain both the shielding of the 
companion and the scattering. Such a model would resemble magnetically driven wind mod- 
els for the broad-line region in AGN (e.g. Emmering et al., 1992, Konigl and Kartje 1994). A 
promising possibility is that the reprocessing region consists of matter 'kicked up' at the disk 
edge by the impact of the mass transferring stream (Meyer-Hofmeister et al. 1997, Armitage 
and Livio 1998, Spruit et al. 1998). This produces qualitatively the right dependence of X-ray 
absorption on orbital phase in ADC sources, and the light curves of the so-called supersoft 
sources. 
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7.3 Transients 

Soft X-ray transients (also called X-ray Novae) are believed to be binaries similar to the other 
LMXB, but somehow the accretion is episodic, with very large outbursts recurring on time 
scales of decades (sometimes years). Most of these transients turn out to be black hole can- 
didates (see Lewin et al. 1995 for a review). As with the Dwarf Novae, the time dependence 
of the accretion can in principle be exploited in transients to derive information on the disk 
viscosity, assuming that the outburst is caused by an instability in the disk. The closest rela- 
tives of soft transients among the White Dwarf plus main sequence star systems are probably 
the WZ Sge stars (van Paradijs and Verbunt 1984, KuuUcers et al. 1996), which show (in the 
optical) similar outbursts with similar recurrence times (cf. Warner 1987, O'Donoghue et al. 
1991). Like the soft transients, they have low mass ratios (q < 0.1). For a given angular mo- 
mentum loss, systems with low mass ratios have low mass transfer rates, so the speculation 
is that the peculiar behavior of these systems is somehow connected in both cases with a low 
mean accretion rate. 



7.3.1 Transients in quiescence 

X-ray transients in quiescence (i.e. after an outburst) usually show a very low X-ray lumi- 
nosity. The mass transfer rate from the secondary in quiescence can be inferred from the 
optical emission. This shows the characteristic 'hot spot', known from other systems to be 
the location where the mass transferring stream impacts on the edge of an accretion disk (e.g. 
van Paradijs and McClintock 1995). These observations thus show that a disk is present in 
quiescence, while the mass transfer rate can be measured from the brightness of the hot spot. 
If this disk were to extend to the neutron star with constant mass flux, the predicted X-ray 
luminosity would be much higher than observed. This has traditionally been interpreted as a 
consequence of the fact that in transient systems, the accretion is not steady. Mass is stored 
in the outer parts and released by a disk instability (e.g. King 1995, Meyer-Hofmeister and 
Meyer 1999) producing the X-ray outburst. During quiescence, the accretion rate onto the 
compact object is much smaller than the mass transfer from the secondary to the disk. 



7.4 Disk Instability 

The most developed model for outbursts is the disk instability model of Osaki (1974), Hoshi 
(1979), Smak (1971, 1984), Meyer and Meyer-Hofmeister (1981), see also King (1995), Os- 
aki (1993). In this model the instability that gives rise to cyclic accretion is due to a tempera- 
ture dependence of the viscous stress. In any local process that causes an effective viscosity, 
the resulting a- parameter will be a function of the main dimensionless parameter of the disk, 
the aspect ratio H/r. If this is a sufficiently rapidly increasing function, such that a is large in 
hot disks and low in cool disks, an instability results by the following mechanism. Suppose 
we start the disk in a stationary state at the mean accretion rate. If this state is perturbed by a 
small temperature increase, a goes up, and by the increased viscous stress the mass flux M 
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increases. By p5\ this increases the disk temperature further, resulting in a runaway to a hot 
state. Since M is larger now than the average, the disk empties partly, reducing the surface 



density and the central temperature (eq. 40 1. A cooling front then transforms the disk to a 
cool state with an accretion rate below the mean. The disk in this model switches back and 
forth between hot and cool states. By adjusting the value of a in the hot and cool states, or by 
adjusting the functional dependence of a on H/r, outbursts are obtained that agree reason- 
ably with the observations of soft transients (Lin and Taam 1984, Mineshige and Wheeler, 
1989). A rather strong increase of a with H/r is needed to get the observed long recurrence 
times. 

Another possible mechanism for instability has been found in 2-D numerical simulations 
of accretion disks (Blaes and Hawley 1988, Rozyczka and Spruit 1993). The outer edge of a 
disk is found, in these simulations, to become dynamically unstable to an oscillation which 
grows into a strong eccentric perturbation (a crescent shaped density enhancement which 
rotates at the local orbital period). Shock waves generated by this perturbation spread mass 
over most of the Roche lobe; at the same time the accretion rate onto the central object is 
strongly enhanced. This process is different from the Smak-Osaki-Hoshi mechanism, since it 
requires 2 dimensions, and does not depend on the viscosity (instead, the internal dynamics 
in this instability generates the effective viscosity that causes a burst of accretion). 



7.5 Other Instabilities 

Instability to heating/cooling of the disk can be the due to several effects. The cooling rate of 
the disk, if it depends on temperature in an appropriate way, can cause a thermal instability 
like that in the interstellar medium. Other instabilities may result from the dependence of 
viscosity on conditions in the disk. For a general treatment see Piran (1978), for a shorter 
discussion see Treves et al. (1988). 



8 Sources of Viscosity 

The high Reynolds number of the flow in accretion disks (of the order lO''* in the outer 
parts of a CV disk) would, to most fluid dynamicists, seem an amply sufficient condition 



for the occurrence of hydrodynamic turbulence (see also the discussion in section 4.7 1. A 
theoretical argument against such turbulence often used in astrophysics (Kippenhahn and 
Thomas 1981,Pringle 1981) is that in cool disks the gas moves almost on Kepler orbits, which 
are quite stable (except for the orbits that get close to the companion or near a black hole). 
This stability is related to the known stabilizing effect that rotation has on hydrodynamical 
turbulence (Bradshaw 1969, Lesur & Longaretti 2005). A (not very strong) observational 
argument is that hydrodynamical turbulence as described above would produce an a that does 
not depend on the nature of the disk, so that all objects should have the same value, which 
is not what observations show. From the modeling of CV outbursts one knows, for example, 
that a probably increases with temperature (more accurately, with H/r, see previous section). 
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Also, there are indications from the inferred life times and sizes of protostellar disks (Strom 
et al. 1993) that a may be rather small there, ^ 10^^, whereas in outbursts of CV's one infers 
values of the order 0.1 — 1. 

Among the processes that have been proposed repeatedly as sources of viscosity is convec- 
tion due to a vertical entropy gradient (e.g. Kley et al. 1993), which may have some limited 
effect in convective parts of disks. Rotation rate varies somewhat with height in a disk, due 
to the radial temperature gradient ('baroclinicity'). Modest amounts of turbulence have been 
reported from instability of this form of differential rotation (Klahr & Bodenheimer 2003, 
Lesur and Papaloizou 2010). 

Another class are waves of various kinds. Their effect can be global, that is, not reducible 
to a local viscous term because by traveling across the disk they can communicate torques 
over large distances. For example, waves set up at the outer edge of the disk by tidal forces can 
travel inward and by dissipating there can effectively transport angular momentum outward 
(e.g. Narayan et al. 1987, Spruit et al. 1987). A nonlinear version of this idea are selfsimilar 
spiral shocks, observed in numerical simulations (Sawada et al. 1987) and studied analytically 
(Spruit 1987). Such shocks can produce accretion at an effective a of 0.01 in hot disks, but 
are probably not very effective in disks as cool as those in CV's and XRB. 

A second non-local mechanism is provided by a magnetically accelerated wind originating 
from the disk surface (Blandford 1976, Bisnovatyi-Kogan and Ruzmaikin 1976, Lovelace 
1976, for an introduction see Spruit 1996). In principle, such winds can take care of all the 
angular momentum loss needed to make accretion possible in the absence of a viscosity 
(Blandford 1976, Konigl 1989). The attraction of this idea is that magnetic winds are a strong 
contender for explaining the strong outflows and jets seen in some protostellar objects and 
AGN. It is not yet clear however if, even in these objects, the wind is actually the main source 
of angular momentum loss. 

In sufficiently cool or massive disks, selfgravitating instabilities of the disk matter can 
produce internal friction. Paczynski (1978) proposed that the resulting heating would limit 
the instability and keep the disk in a well defined moderately unstable state. The angular 
momentum transport in such a disk has been studied numerically (e.g. Gammie 1997, Ostriker 
et al. 1999). Disks in CVs and XRB are too hot for selfgravity to play a role, but it can be 
important in protostellar disks (cf. Rafikov 2009). 



8.1 magnetic viscosity 

Magnetic forces can be very effective at transporting angular momentum. If it can be shown 
that the shear flow in the disk produces some kind of small scale fast dynamo process, that is, 
some form of magnetic turbulence, an effective a ^ 0(1) would be expected (Shakura and 
Sunyaev 1973). Numerical simulations of initially weak magnetic fields in accretion disks 
show that this does indeed happen in sufficiently ionized disks (Hawley et al. 1995, Bran- 
denburg et al. 1995, Balbus 2003). These show a small scale magnetic field with azimuthal 
component dominating (due to stretching by differential rotation). The angular momentum 
transport is due to magnetic stresses, the fluid motions induced by the magnetic forces con- 
tribute only little to the angular momentum transport. In a perfectly conducting plasma this 
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turbulence can develop from an arbitrarily small initial field through magnetic shear insta- 
bility (also called magnetorotational instability, Velikhov 1959, Chandrasekhar 1961, Balbus 
and Hawley 1991). The significance of this instability is that it shows that at large conductiv- 
ity accretion disks must be magnetic. 

The actual form of the highly time dependent small scale magnetic field which develops 
can only be found from numerical simulations. The simplest case thought to be representative 
of the process considers a vertically unstratified disk (component of gravity perpendicular 
to the disk ignored), of which only a radial extent of a few times H is included (cf. the 
causality argument above). Different simulations have yielded somewhat different values for 
the effective viscosity, contrary to the expectation that at least in this simplest form the process 
should yield a unique value. 

The nature of these differences has been appreciated only recently. Fromang et al. (2007a) 
show that the results do not converge with increasing numerical resolution (the effective a 
found appears to decrease indefinitely). The result also depends on the magnetic Prandtl 
number P^ — v/rjni, the ratio of viscosity to magnetic diffusivity Tj^. No turbulence is found 
for Pm < 1 Fromang et al. (2007b). The significance of these findings is not quite obvious 
yet, but they clearly contradict the commonly assumed 'cascade' picture of MHD turbulence 
(taken over from what is assumed in hydrodynamic turbulence, where the large scale behavior 
of turbulence appears to converge with numerical resolution). 

Perhaps the behavior found in this case is related to the highly symmetric nature of the 
idealized problem. This is suggested by the finding (Davis et al. 2010, Shi et al. 2010) that 
effective viscosity appears to converge again with increasing numerical resolution when ver- 
tical stratification is included. By inference, the process defining the state of magnetic turbu- 
lence is different in this case from that in the unstratified case. It can plausibly be attributed 
to magnetic buoyancy instabilities (Shi et al. 2010). This would make the process somewhat 
analogous to the mechanism operating the solar cycle (Spruit 2010). Magnetic turbulence 
also behaves different if the net vertical magnetic flux crossing the simulated box does not 
vanish. (This flux is a conserved quantity set by the initial conditions). Magnetic turbulence 
then develops more easily (at lower numerical resolution), and appears to converge as the 
numerical resolution is increased. 

The consequences of these new findings for angular momentum transport in disks are still 
to be settled (cf the insightful discussion in Lesur and Ogilvie 2008). 



8.2 viscosity in radiatively supported disks 

A disk in which the radiation pressure fr dominates must be optically thick (otherwise the 
radiation would escape). The radiation pressure then adds to the total pressure. The pressure 
is larger than it would be, for a given temperature, if only the gas pressure were effective. If 



the viscosity is then parametrized by ( 19 1, it turns out (Lightman and Eardley, 1974) that the 
disk is locally unstable. An increase in temperature increases the radiation pressure, which 
increases the viscous dissipation and the temperature, leading to a runaway. This has raised 
the question whether the radiation pressure should be included in the sound speed that enters 



expression ( 19 1. If it is left out, a lower viscosity results, and there is no thermal-viscous run- 
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away. Without knowledge of the process causing the effective viscous stress, this question can 
not be answered. Sakimoto and Coroniti (1989) have shown, however, that if the stress is due 
to some form of magnetic turbulence, it most likely scales with the gas pressure alone, rather 
than the total pressure. Now that it seems likely, from the numerical simulations, that the 
stress is indeed magnetic, there is reason to believe that in the radiation pressure-dominated 
case the effective viscosity will scale as v ^ aPg/{pQ), making such disks thermally stable, 
as indicated by the results of Hirose et al. (2009). 



9 Beyond thin disks 

Ultimately, much of the progress in developing useful models of accretion disks will depend 
on detailed numerical simulations in 2 or 3 dimensions. In the disks one is interested in, there 
is usually a large range in length scales (in LMXB disks, from less than the 10 km neutron star 
radius to the more than 10^ km orbital scale). Correspondingly, there is a large range in time 
scales that have to be followed. This not technically possible at present and in the foreseeable 
future. In numerical simulations one is therefore limited to studying in an approximate way 
aspects that are either local or of limited dynamic range in r, t (for examples, see Hawley 
1991, Armitage 1998, De Villiers et al. 2003, Hirose et al. 2004). For this reason, approaches 
have been used that relax the strict thin disk framework somewhat without resorting to full 
simulations. Some of the physics of thick disks can be included in a fairly consistent way in 
the 'slim disk' approximation (Abramowicz et al., 1988). The so-called Advection Dominated 
Accretion Flows (ADAFs) are related to this approach (for a review see Yi 1998). They are 
discussed in sections [TT] [T2] [T3] below). 



9.1 Boundary layers 

In order to accrete onto a star rotating at the rate i2*, the disk matter must dissipate an amount 

of energy given by 

GMM , , , ,,2 

-^^[\-a,/£lk{R)f. (59) 

The factor in brackets measures the kinetic energy of the matter at the inner edge of the disk 
(r = R), in the frame of the stellar surface. Due to this dissipation the disk inflates into a 'belt' 
at the equator of the star, of thickness D and radial extent of the same order. Equating the 



radiation emitted from the surface of this belt to ( 59 1 one gets for the surface temperature Tjb 



of the belt, assuming optically thick conditions and a slowly rotating star (Q^,/Qk <C 1): 

GMM A 

Y- = (^vTI (60) 



To find the temperature inside the belt and its thickness, use eq. (39 1. The value of the surface 
temperature is higher, by a factor of the order (R/DY'^, than the simplest thin disk estimate 
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([35] ignoring the factor /). In practice, this works out to a factor of a few. The surface of 
the belt is therefore not very hot. The situation is quite different if the boundary layer is not 
optically thick (Pringle and Savonije 1979). It then heats up to much higher temperatures. 

Analytical methods to obtain the boundary layer structure have been used by Regev and 
Hougerat (1988), numerical solutions of the slim disk type by Narayan and Popham (1993), 
Popham (1997), 2-D numerical simulations by Kley (1991). These considerations are primar- 
ily relevant for CV disks; in accreting neutron stars, the dominant effects of radiation pressure 
have to be included. More analytic progress on the structure of the boundary layer between a 
disk and a neutron star and the way in which it spreads over the surface of the star has been 
reported by Inogamov and Sunyaev (1999). 



10 Radiative efficiency of accretion disks 

In a thin accretion disk, the time available for the accreting gas to radiate away the energy 
released by the viscous stress is the accretion time, 

' 0\ (61) 



'ace ~' r^ \ TT 

uQk H 

where a is the dimensionless viscosity parameter, Q}^ the local Keplerian rotation rate, r the 
distance from the central mass, and H the disk thickness (see Frank et al. 2002 or section[6|. 
For a thin disk, H/r <C 1, this time is much longer than the thermal time scale ft ~ l/{aQ) 
(section|6]l. There is then enough time for a local balance to hold between viscous dissipation 
and radiative cooling. For the accretion rates implied in observed systems the disk is then 
rather cool, which then justifies the starting assumption H/r ^ 1. 

This argument is somewhat circular, of course, since the accretion time is long enough for 
effective cooling only if the disk is assumed to be thin to begin with. Other forms of accretion 
disks may exist, even at the same accretion rates, in which the cooling is ineffective compared 
with that of standard (geometrically thin, optically thick) disks. In the following sections, we 
consider such forms of accretion and the conditions under which are to be expected. 

Since radiatively inefficient disks tend to be thick, H/r ^ 0{l), they are sometimes called 
'quasi-spherical'. However, this does not mean that a spherically symmetric accretion model 
would be a reasonable approximation. The crucial difference is that the flow has angular 
momentum. The inward flow speed is governed by the rate at which angular momentum can 
be transferred outwards, rather than by gravity and pressure gradient as in the Bondi accretion 
problem mentioned above. With H/r ^ ^{l), the accretion time scale, face ^ l/{aQ) is still 
longer than the accretion time scale in the spherical case, l/Q, (unless the viscosity parameter 
a is as large as 0(1)). The dominant velocity component is azimuthal rather than radial, and 
the density and optical depth are much larger than in the spherical case. 

It turns out that there are two kinds of radiatively inefficient disks, the optically thin and 
optically thick varieties. A second distinction occurs because accretion flows are different for 
central objects with a solid surface (neutron stars, white dwarfs, main sequence stars, planets), 
and those without (i.e. black holes). Start with optically thick flows. 
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11 Radiation supported radiatively inefficient accretion 

If the energy loss by radiation is small, the gravitational energy release Wgrav ~ GM/{2r) is 
converted into enthalpy of the gas and radiation fielqj 

''?^.i[JLp^ + 4PJ. ,62, 

2 r p y- 1 " 

where an ideal gas of constant ratio of specific heats y has been assumed, and P^ ~ \aj'' 
is the radiation pressure. In terms of the virial temperature Tm„ — GMjijXr), and assuming 



7 = 5/3, appropriate for a fully ionized gas (see section 2.1 1, this can be written as 



^ = [5 + 8^^]-^ (63) 

Thus, for radiation pressure dominated accretion, P^ ^ /g, the temperature is much less than 
the virial temperature: much of the accretion energy goes into photon production instead of 
heating. By hydrostatic equilibrium the disk thickness is given by (cf section |4~2| 

i/«[(Pg+P,)/p]'/V^, (64) 



With (|63]l this yields 

H/r^O{\). (65) 

In the limit P^ ^ Pg, the flow is therefore geometrically thick. Radiation pressure then supplies 
a non-negligible fraction of the support of the gas in the radial direction against gravity (the 
remainder being provided by rotation). 



For Pr > ^g, (|62|i yields 

GM 4 nT^ 

(66) 
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(67) 

Hence 

(68) 

where Fe = LE/{4-7ir^) is the local Eddington flux. Since H/r w 1, a radiatively inefficient, 
radiation pressure dominated accretion flow has a luminosity of the order of the Eddington 
luminosity. 

The temperature depends on the accretion rate and the viscosity v assumed. The accretion 



rate is of the order M ^ 37i:vZ (cf. eq. (29i), where L — J pdzis the surface mass density. In 



units of the Eddington rate Me, eq. ([8jl, we get 



' This assumes that a fraction ~ 0.5 of the gravitational potential energy stays in the flow as orbital kinetic 
energy. This is only an approximation, see also section! 121 
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in = M/Me « vp k/c, (69) 

where H/r sa 1 has been usecrl Assume that the viscosity scales with the gas pressure: 

v = a-^, (70) 

instead of the total pressure Pi +Pg. This is the form that is likely to hold if the angular 
momentum transport is due to a small-scale magnetic field (Sakimoto and Coroniti, 1989, 



Turner 2004). Then with (66 1, (69 1 we have (up to a numerical factor of ^(1)) 



,^^ iGM)y^ mc 

or 

r«10**r-'/'(r/rg)3/iOmi/5^ (72) 

where r — lO^rg and rg — 2GM/c^ is the gravitational radius of the accreting object, and 
the electron scattering opacity of 0.3 cm^/g has been assumed. The temperatures expected 
in radiation supported advection dominated flows are therefore quite low compared with the 
virial temperature [If the viscosity is assumed to scale with the total pressure instead of Pg, the 
temperature is even lower] . The effect of electron-positron pairs can be neglected (Schultz and 
Price, 1985), since they are present only at temperatures approaching the electron rest mass 
energy, T > lO^K. 

In order for the flow to be radiation pressure and advection dominated, the optical depth 
has to be sufficiently large so the radiation does not leak out. The energy density in the flow, 
vertically integrated is of the order 

E ss aT^H, (73) 



and the energy loss rate per cm of disk surface is given by ( 68 1. The cooling time is therefore. 



fc = E/F = 3t///c. (74) 

This is to be compared with the accretion time, which can be written in terms of the mass in 
the disk at radius r, of the order 2nr^E, and the accretion time: 

face ~ 2nr^E/M. (75) 

This yields 

K 4 R 

fcAacc « M = -m-, (76) 

Ttrc rj r 

(where a factor Ji/IH /r ~ 0(1) has been neglected). Since r > R, this shows that accretion 
has to be around the Eddington rate or larger in order to be both radiation- and advection- 
dominated. 



The definition of Me differs by factors of order unity between different authors. It depends on the assumed 
efficiency 77 of conversion of gravitational energy GM/R into radiation. For accretion onto black holes a more 
realistic value is of order 7] = 0.1, for accretion onto neutron stars 77 r; 0.4, depending on the radius of the 
star. 
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This condition can also be expressed in terms of the so-called trapping radius r^ (e.g. Rees 
1978). Equating face and ^c yields 

r^/R sa 4m. (77) 

Inside r^, the flow is advection dominated: the radiation field produced by viscous dissipation 
stays trapped inside the flow, instead of being radiated from the disk as happens in a stan- 
dard thin disk. Outside the trapping radius, the radiation field cannot be sufficiently strong to 
maintain a disk with H /r ~ 1, it must be a thin form of disk instead. Such a thin disk can still 
be radiation-supported (i.e. Pr ^ Pg), but it can not be advection dominated. 

Flows of this kind are called 'radiation supported accretion tori' (or radiation tori, for 
short) by Rees et al. 1982. They must accrete at a rate above the Eddington value to exist. 
The converse is not quite true: a flow accreting above Eddington is an advection dominated 
flow, but it need not necessarily be radiation dominated. Advection dominated optically thick 



accretion flows exist in which radiation does not play a major role (see section 12.1 1. 

That an accretion flow above Me is advection dominated, not a thin disk, also follows from 
the fact that in a thin disk the energy dissipated must be radiated away locally. Since the local 
radiative flux can not exceed the Eddington energy flux fg, the mass accretion rate in a thin 
disk can not significantly exceed the Eddington value (|8]l. 

The gravitational energy, dissipated by viscous stress in differential rotation and advected 
with the flow, ends up on the central object. If this is a black hole, the photons, particles 
and their thermal energy are conveniently swallowed at the horizon, and do not react back 
on the flow. Radiation tori are therefore mostly relevant for accretion onto black holes. They 
are convectively unstable (Bisnovatyi-Kogan and Blinnikov 1977): the way in which energy 
is dissipated, in the standard a -prescription, is such that the entropy (entropy of radiation. 



r^ T^ Ip) decreases with height in the disk. Numerical simulations (see section 14 1 show the 
effects of this convection. 



11.1 Super-Eddington accretion onto black holes 

As the accretion rate onto a black hole is increased above Me, the trapping radius moves 
out. The total luminosity increases only slowly, and remains of the order of the Eddington 
luminosity. Such supercritical accretion has been considered by Begelman and Meier (1982, 
see also Wang and Zhou 1999); they show that the flow has a radially self-similar structure. 

Abramowicz et al. (1988, 1989) studied accretion onto black holes at rates near Mg. They 
used a vertically-integrated approximation for the disk, but included the advection terms. 
The resulting models were called 'slim disks'. They show how with increasing accretion 
rate, a standard thin Shakura-Sunyaev disk turns into a radiation-supported advection flow. 
The nature of the transition depends on the viscosity prescription used, and can show a non- 
monotonic dependence of M on surface density T, (Honma et al. 1991). This suggests the 
possibility of instability and cyclic behavior of the inner disk near a black hole, at accretion 
rates near and above Me (for an application to GRS 1915H-105 see Nayakshin et al., 1999). 
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11.2 Super-Eddington accretion onto neutron stars 

In the case of accretion onto a neutron star, the energy trapped in the flow, plus the remaining 
orbital energy, settles onto its surface. If the accretion rate is below Me, the energy can be 
radiated away by the surface, and steady accretion is possible. A secondary star providing the 
mass may, under some circumstances, transfer more than Me, since it does not know about 
the neutron star's Eddington value. The outcome of this case is still somewhat uncertain; it 
is generally believed on intuitive grounds that the 'surplus' (the amount above Me) somehow 
gets expelled from the system. 

One possibility is that, as the transfer rate is increased, the accreting hot gas forms an 
extended atmosphere around the neutron star like the envelope of a giant. If it is large enough, 
the outer parts of this envelope are partially ionized. The opacity in these layers, due to atomic 
transitions of the CNO and heavier elements, is then much higher than the electron scattering 
opacity. The Eddington luminosity based on the local value of the opacity is then smaller 
than it is near the neutron star surface. Once an extended atmosphere forms, the accretion 
luminosity is thus large enough to drive a wind from the envelope (see Kato 1997, where the 
importance of this effect is demonstrated in the context of Novae). 

This scenario is somewhat dubious however, since it assumes that the mass transferred 
from the secondary continues to reach the neutron star and generate a high luminosity there. 
This is not at all obvious, since the mass transferring stream may instead dissipate inside the 
growing envelope of the neutron star. The result of this could be a giant (more precisely, a 
Thorne-Zytkow star), with a steadily increasing envelope mass. Such an envelope is likely 
to be large enough to engulf the entire binary system, which then develops into a common- 
envelope (CE) system. The envelope mass is then expected to be ejected by CE hydrodynam- 
ics (for reviews see Taam 1994, Taam and Sandquist 2000). 

A more speculative proposal, suggested by the properties of SS 433, is that the 'surplus 
mass' is ejected in the form of jets. The binary parameters of Cyg X-2 are observational 
evidence for mass ejection in super-Eddington mass transfer phases (King and Ritter 1999, 
Rappaport and Podsiadlowski 2000, King and Begelman 1999). 



12 ADAF Hydrodynamics 

The hydrodynamics of radiatively inefficient flows (or 'advection dominated accretion flows') 
can be studied by starting, at a very simple level, with a generalization of the thin disk equa- 
tions. Making the assumption that quantities integrated over the height z of the disk give a 
fair representation (though this is justifiable only for thin disks), and assuming axisymmetry, 
the problem reduces to a one-dimensional time-dependent one. Further simplifying this by 
restriction to a steady flow yields the equations 

InrZvr =M = est, (78) 

rEvrdriQr^) = drivEr^drQ), (79) 
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(80) 



EvrTd,.S = q+-q-, (81) 

where S is the specific entropy of the gas, Q the local rotation rate, now different from the 
Keplerian rate XSk, and 

q+= [g^dz q-=U\vFAz (82) 



are the height-integrated viscous dissipation rate and radiative loss rate, respectively. In the 
case of thin disks, equations ( |78]l an d (79 1 are unchanged, but (80i simplifies to i2^ = i2^, i.e. 
the rotation is Keplerian, while ([8 1} simplifies to q^ = q^, expressing local balance between 
viscous dissipation and cooling. The left hand side of ( [8T| l describes the radial advection of 
heat, and is perhaps the most important deviation from the thin disk equations at this level 
of approximation (hence the name advection dominated flows). The characteristic properties 
are seen most clearly when radiative loss is neglected altogether, q^ = 0. The equations are 
supplemented with expressions for v and q^: 



V = acl/Q.^; q+ = {rdrQfvE. 



(83) 



If a is taken constant, q =0, and an ideal gas is assumed with constant ratio of specific 
heats, so that the entropy is given by 



S — Cvln{p/p'> 



(84) 



then equations (78 i-(|8T|l have no explicit length scale in them. This means that a special so- 



called self-similar solution exists, in which all quantities are powers of r. Such self-similar 
solutions have apparently been described first by Gilham (1981), but re-invented several times 
(Spruit et al. 1987; Narayan and Yi, 1994). The dependences on r are 



Q 



-3/2. 



H' 



' r 
.-1 



3/2 



In the limit a ^ 1, one finds 
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The precise from of these expressions depends somewhat on the way in which vertical inte- 



grations such as in ( 82 1 are done (which are only approximate). 



The self-similar solution can be compared with numerical solutions of eqs. ( 78 1-( 8 1 1 with 
appropriate conditions applied at inner (r;) and outer (ro) boundaries (Nakamura et al. 1996, 
Narayan et al. 1997). The results show that the self-similar solution is valid in an intermediate 
regime n ^ r ^ Tq. That is, the solutions of ([78]l-([8T]l approach the self-similar solution far 
from the boundaries, as is characteristic of self-similar solutions. 

The solution exists only if 1 < 7 < 5/3, a condition satisfied by all ideal gases. As 7I 1, 
the disk temperature and thickness vanish. This is understandable, since a 7 close to 1 means 
that the particles making up the gas have a large number of internal degrees of freedom. In 
thermal equilibrium the accretion energy is shared between all degrees of freedom, so that 
for a low 7 less is available for the kinetic energy (temperature) of the particles. 

Second, the rotation rate vanishes for 7—)- 5/3. As in the case of spherical accretion no 
accreting solutions exist for 7 > 5/3 (cf. sectionl2]). Since a fully ionized gas has 7=5/3, it is 
a relevant value for the hot, ion supported accretion flows discussed below. Apparently, steady 
advection dominated accretion can not have angular momentum in this case. To see how 
comes about consider the entropy of the flow. For accretion to take place in a rotating flow, 
there has to be friction and increase of entropy. Accretion is then necessarily accompanied 



by an inward increase of entropy. In a radially self similar flow, the scalings ( 85 1, ( 86 1 yield 
P ^ r^^'^, so entropy scales as 5 ^ InP/p^ ^ (5 — 37)lnr. This increases with decreasing 
distance r only for 7 < 5/3, and an accretion flow with 7 = 5/3 cannot be both rotating and 
adiabatic. (With energy losses by radiation the constraint on 7 disappears again). 

The question then arises how an adiabatic flow with 7 = 5/3 will behave if one starts it as 
a rotating torus around a black hole. In the literature, this problem has been circumvented by 
arguing that real flows would have magnetic fields in them, which would change the effective 
compressibility of the gas. Even if a magnetic field of sufficient strength is present, however, 
(energy density comparable to the gas pressure) the effective 7 is not automatically lowered. 
If the field is compressed mainly perpendicular to the field lines, for example, the effective 
7 is closer to 2. Also, this does not solve the conceptual problem what would happen to a 
rotating accretion flow consisting of a more weakly magnetized ionized gas. 



12.0.1 The case of the vanishing rotation rate 

This conceptual problem has been solved by Ogilvie (1999), who showed how a gas cloud 
initially rotating around a point mass settles to the slowly rotating self-similar solutions of the 
steady problem discussed above. He constructed similarity solutions to the time dependent 



version of eqs (78 1-(81 1, in which distance and time occur in the combination r/f^". This 



solution describes the asymptotic behavior (in time) of a viscously spreading disk, analogous 



to the viscous spreading of thin disks (see section 4.3 1. As in the thin disk case, all the mass 
accretes asymptotically onto the central mass, while all the angular momentum travels to 
infinity together with a vanishing amount of mass. For all 7 < 5/3, the rotation rate at a 
fixed r tends to a finite value as f — >^ 00, but for 7=5/3 it tends to zero. The size of the 
slowly -rotating region expands as r ^ f^' ^.The typical slow rotation of ADAFs at 7 near 5/3 
is thus a real physical property. In such a flow the angular momentum gets expelled from the 
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inner regions almost completely, and the accretion flow becomes purely radial, as in Bondi 
accretion. 



12.1 Other optically thick accretion flows 



The radiation-dominated flows discussed in section 11 are not the only possible optically 
thick advection dominated flows. From the discussion of the hydrodynamics, it is clear that 
disk-like (i.e. rotating) accretion is possible whenever the ratio of specific heats is less than 
5/3. A radiation supported flow satisfies this requirement since radiation pressure scales with 
volume like a gas with y = 4/3, but it can also happen in the absence of radiation if en- 
ergy is taken up in the gas by internal degrees of freedom of the particles. Examples are the 
rotational and vibrational degrees of freedom in molecules, and the energy associated with 
dissociation and ionization. If the accreting object has a gravitational potential not too far 
from the 2.3 + 13.6 eV per proton for dissociation plus ionization, a gas initially consisting 
of molecular hydrogen can stay bound at arbitrary accretion rates. This translates into a limit 
M/MqRq/R < 0.01. This is satisfied approximately by the giant planets, which are believed 
to have gone through a phase of rapid adiabatic gas accretion (e.g. Podolak et al. 1993). 

A more remotely related example is the core-collapse supernova. The accretion energy of 
the core mass falling onto the growing proto-neutron star at its center is lost mostly to internal 
degrees of freedom represented by photodisintegration of the nuclei. If the pre-collapse core 
rotates sufficiently rapidly, the collapse will form an accretion torus (inside the supernova 
envelope), with properties similar to advection dominated accretion flows (but at extreme 
densities and accretion rates, by X-ray binary standards). Such objects have been invoked as 
sources of Gamma-ray bursts (Woosley 1993, Paczyriski 1998, Popham et al. 1999). 

A final possibility for optically thick accretion is through neutrino losses. If the tempera- 
ture and density near an accreting neutron star become large enough, additional cooling takes 
place through neutrinos (as in the cores of giants). This is relevant for the physics of Thome- 
Zytkow stars (neutron stars or black holes in massive supergiant envelopes, cf. Bisnovatyi- 
Kogan and Lamzin 1984, Cannon et al. 1992), and perhaps for the spiral-in of neutron stars 
into giants (Chevalier 1993, see however Taam & Sandquist 2000). 



13 Optically thin radiatively inefficient flows (ISAFs) 

The optically thin case has received most attention, because of the promise it holds for ex- 
plaining the (radio to X-ray) spectra of X-ray binaries and the central black holes in galaxies, 
including our own. For a review see Yi (1999). This kind of flow occurs if the gas is optically 
thin, and radiation processes are sufficiently weak. The gas then heats up to near the virial 
temperature. Near the last stable orbit of a black hole, this is of the order 100 MeV, or 10^^ 
K. At such temperatures, a gas in thermal equilibrium would radiate at a fantastic rate, even 
if it were optically thin, because the interaction between electrons and photons becomes very 
strong already near the electron rest mass of 0.5 MeV. In a remarkable early paper, Shapiro, 
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Lightman and Eardley (1976) noted that this, however, is not what will happen in an optically 
thin plasma accreting on a hole. They showed that instead thermal equilibrium between ions 
and electrons breaks down and a two-temperature plasma forms. We call such a flow an ion 
supported accretion flow (ISAF), following the nomenclature suggested by Rees et al. (1982). 
The argument is as follows. 

Suppose that the energy released by viscous dissipation is distributed equally among the 
carriers of mass, i.e. mostly to the ions and ^ 1/2000 to the electrons. Most of the energy 
then resides in the ions, which radiate very inefficiently (their high mass prevents the rapid 
accelerations that are needed to produce electromagnetic radiation). Their energy is transfered 
to the electrons by Coulomb interactions. These interactions are slow, however, under the 
conditions mentioned. They are slow because of the low density (on account of the assumed 
low optical thickness), and because they decrease with increasing temperature. The electric 
forces that transfer energy from an ion to an electron act only as long as the ion is within 
the electron's Debye sphere (e.g. Spitzer, 1965). The interaction time between proton and 

— 1/2 

electron, and thus the momentum transfered, therefore decrease as 1/vp ^ Tp where Tp is 
the proton temperature. 

In this way, an optically thin plasma near a compact object can be in a two-temperature 
state, with the ions being near the virial temperature and the electrons, which are doing the 
radiating, at a much lower temperature around 50-200 keV. The energy transfer from the 
gravitational field to the ions is fast (by some form of viscous or magnetic dissipation, say), 
from the ions to the electrons slow, and finally the energy losses of the electrons fast (by 
synchrotron radiation in a magnetic field or by inverse Compton scattering off soft photons). 
Such a flow would be radiatively inefficient since the receivers of the accretion energy, the 
ions, get swallowed by the hole before having a chance to transfer their energy to the elec- 
trons. Most of the accretion energy thus gets lost into the hole, and the radiative efficiency is 
much less than for a cool disk. The first disk models which take into account this physics of 
advection and a two-temperature plasma were developed by Ichimaru (1977). 

It is clear from this description that both the physics of such flows and the radiation spec- 
trum to be expected depend critically on the details of the ion-electron interaction and radi- 
ation processes assumed. This is unlike the case of the optically thick advection dominated 
flows, where gas and radiation are in approximate thermodynamic equilibrium. This is a 
source of uncertainty in the application of ISAFs to observed systems, since their radiative 
properties depend on poorly known quantities such as the strength of the magnetic field in 
the flow. 

The various branches of optically thin and thick accretion flows are summarized in figure 
HI Each defines a relation between surface density Z (or optical depth T = kT.) and accretion 
rate. ISAFs require low densities, which can result either because of low accretion rates, or 
large values of the viscosity parameter. The condition that the cooling time of the ions by 
energy transfer to the electrons is longer than the accretion time yields a maximum accretion 
rate (Rees et al. 1982), 

m<,a^, (91) 

where m is the accretion rate in units of the Eddington value. If a w 0.05 as suggested by 
current simulations of magnetic turbulence, the maximum accretion rate would be a few 10^^. 
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Fig. 4 Branches of advection-dominated and thin disks for two values or the viscosity parameter a, as func- 
tions of accretion rate and (vertical) optical depth of the flow (schematic, after Chen et al. 1995, Zdziarski 
1998). Optically thin branches are the ISAF and SLE (Shapiro-Lightman-Eardley) solutions, optically thick 
ones the radiation dominated ('slim disk' or 'radiation torus') and SS (Shakura-Sunyaev or standard thin 
disk). Advection dominated are the ISAF and the radiation torus, geometrically thin are the SLE and SS. The 
SLE solution is a thermally unstable branch. 



If ISAFs are to be applicable to systems with higher accretion rates, such as Cyg X-1 for 
example, the viscosity parameter must be larger, on the order of 0.3. 



13.1 Application: hard spectra in X-ray binaries 



In the hard state, the X-ray spectrum of black hole and neutron star accreters is characterized 
by a peak in the energy distribution (v/v or EF{E)) at photon energies around 100 keV. This 
is to be compared with the typical photon energy of ^ 1 keV expected from a standard opti- 
cally thick thin disk accreting near the Eddington limit. The standard, and by far most likely 
explanation is that the observed hard photons are softer photons (around 1 keV) that have 
been up-scattered through inverse Compton scattering on hot electrons. Fits of such Comp- 
tonized spectra (e.g. Sunyaev and Titarchuk 1980, Zdziarski 1998 and references therein) 
yield an electron scattering optical depth around unity and an electron temperature of 50-100 
keV. The scatter in these parameters is rather small between different sources. The reason 
may lie in part in the physics of Comptonization, but this is not the only reason. Something 
in the physics of the accretion flow keeps the Comptonization parameters fairly constant as 
long as it is in the hard state. ISAFs have been applied with some success in interpreting 
XRB. They can produce reasonable X-ray spectra, and have been used in interpretations of 
the spectral-state transitions in sources like Cyg X-1 (Esin et al. 1998 and references therein). 
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An alternative to the ISAF model for the hard state in sources like Cyg X-1 and the black 
hole X-ray transients is the 'corona' model. A hot corona (Bisnovatyi-Kogan and Blinnikov 
1976), heated perhaps by magnetic fields as in the case of the Sun (Galeev et al. 1979) could 
be the medium that Comptonizes soft photons radiated from the cool disk underneath. The 
energy balance in such a model produces a Comptonized spectrum within the observed range 
(Haardt and Maraschi 1993). This model has received further momentum, especially as a 
model for AGN, with the discovery of broadened X-ray lines interpreted as indicative of the 
presence of a cool disk close to the last stable orbit around a black hole (Fabian et al. 2002 
and references therein). The very rapid X-ray variability seen in some of these sources is 
interpreted as due magnetic flaring in the corona, like in the solar corona (e.g. Di Matteo et 
al. 1999a). 



13.2 Transition from cool disk to ISAF: truncated disks 

One of the difficulties in applying ISAFs to specific observed systems is the transition from 
a standard geometrically thin, optically thick disk, which must be the mode of mass transfer 
at large distances, to an ISAF at closer range. This is shown by figure p] which illustrates 
the situation at some distance close to the central object. The standard disk and the optically 
thin branches are separated from each other for all values of the viscosity parameter This 
separation of the optically thin solutions also holds at larger distances. Thus, there is no 
plausible continuous path from one to the other, and the transition between the two must be 
due to additional physics that is not included in diagrams like figure l4] 

Circumstantial observational evidence points to the existence of such a transition. The dis- 
tance from the hole where it is assumed to take place is then called the truncation radius. 
The extensive datasets from the black hole candidate Cyg X-1 obtained with the Rossi X-ray 
Timing Explorer (RXTE) have played an important role in the development of the truncated 
disk model. The X-ray spectrum of Cyg X-1 varies (on time scales of days to years) between 
softer and harder states, and the characteristic time scale of its fast variability (milliseconds 
to minutes) correlates closely with these changes. Since all disk time scales decrease with 
distance from the hole, the fast variability is interpreted as due to a process (still to be iden- 
tified in detail) that depends on the size of the truncation radius. Variation in time of this 
radius is then assumed to cause the observed changes. Characteristics of the spectrum that 
find a natural place in this picture are: the slope of the hard X-ray spectrum, the amplitude of 
the so-called 'Compton reflection hump', and the behavior of the Fe K^ fluorescence line at 
6.7 keV; each of these correlates with the characteristic fast variability time scale in an inter- 
pretable (though still somewhat model-dependent) way (Gilfanov et al. 1999, Revnivtsev et 
al. 1999). Very similar correlations have been found in X-ray observations of AGN (Zdziarski 
et al. 1999). 

A promising possibility is that the transition takes place through evaporation. Two dis- 
tinct mechanisms have been elaborated for such evaporation. In the first (Meyer and Meyer- 
Hofmeister 1994, Liu et al. 2002), the evaporation starts at a relatively large distance from the 
hole, where the virial temperature is of the order of 10*' — 10^ K. As in the solar corona, the 
strong decrease of radiative efficiency of gas with temperature in this range produces a hot 
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optically thin corona in contact with the cool disk below, and exchange of mass can take place 
through evaporation and condensation, and the process is mediated by electron heat conduc- 
tion. In this scenario, a corona flow at ^ 10^ K at a distance of several hundred Schwarzschild 
radii transforms into a two-temperature ISAF further in. 



13.2.1 Ion illumination 



evaporation / 



hot protons 



BH 
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Fig. 5 interaction between ion supported accretion flow(ISAF) and cool disk in ttie truncated disk model for 
hard X-ray states in X-ray binaries. Accreting black hole on the left, its hot ions illuminating the cool disk 
create a warm (~ 100 keV surface layer producing the hard X-rays, extending inward and heating up to virial 
temperature by thermal instability ('evaporation') once cooling by soft photons becomes inefficient 



Observations indicate that cool disks can also coexist with a hot, hard X-ray producing 
plasma quite close to the hole (for references see DuUemond and Spruit 2004). At these close 
distances, evaporation must behave differently from the coronal evaporation model, since 
the interaction of a two-temperature plasma with a cool disk is very different from that of 
a plasma at coronal temperatures (Spruit 1997). The energy in an ISAF is in the ions, and 
the electron conduction of heat that drives coronal evaporation unimportant. Moreover, the 
ions penetrate a substantial distance into the cool disk, and the energy they dump is radiated 
away long before they can heat up the disk to virial temperatures. Nevertheless, evaporation 
can still take place in this case, since it turns out that the interaction of the ions with the cool 
disk produces a layer of intermediate temperature (around 100 keV) that becomes thermally 
unstable in the presence of viscous dissipation, and heats up ISAF temperatures (Deufel et 
al. 2000, 2001, 2002, see also Spruit and Deufel 2002). The sequence of events is illustrated 
in figure l5] This model explains both the hard spectra of typical black hole accreters and 
the coexistence of cool and hot plasma that is indicated by the observations (Dullemond and 
Spruit 2004, D'Angelo et al. 2008). 
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13.3 Quiescent galactic nuclei 

For very low accretion rates, such as inferred for the black hole in the center of our galaxy 
(identified with the radio source Sgr A*), the broad band spectral energy distribution of an 
ISAF is predicted to have two humps (Narayan et al. 1995, Quataert et al. 1999). In the X-ray 
range, the emission is due to bremsstrahlung. In the radio range, the flow emits synchrotron 
radiation, provided that the magnetic field in the flow has an energy density order of the gas 
pressure ('equipartition'). Synthetic ISAF spectra can be fitted to the observed radio and X- 
ray emission from Sgr A*. In other galaxies where a massive central black hole is inferred, 
and the center is populated by an X-ray emitting gas of known density, ISAFs would also 
be natural, and might explain why the observed luminosities are so low compared with the 
accretion rate expected for a hole embedded in a gas of the measured density. 



13.4 ISAF-disk interaction: Lithium 

One of the strong predictions of ISAF models, whether for black holes or neutron stars, is 
that the accreting plasma in the inner regions has an ion temperature of 10-100 MeV. Nearby 
is a cool and dense accretion disk feeding this plasma. If only a small fraction of the hot 
ion plasma gets in contact with the disk, the intense irradiation by ions will produce nuclear 
reactions (Aharonian and Sunyaev 1984, Martin et al. 1992). The main effects would be spal- 
lation of CNO elements into Li and Be, and the release of neutrons by various reactions. In 
this context, it is intriguing that the secondaries of both neutron star and black hole accreters 
have high overabundances of Li compared with other stars of their spectral types (Martin et 
al. 1992, 1994a). If a fraction of the disk material is carried to the secondary by a disk wind, 
the observed Li abundances may be accounted for (Martin et al. 1994b). 



14 Outflows 

The energy density in an advection dominated accretion flow is of the same order as the 
gravitational binding energy density GM/r, since a significant fraction of that energy went 
into internal energy of the gas by viscous dissipation, and little of it got lost by radiation. 
The gas is therefore only marginally bound in the gravitational potential. This suggests that 
perhaps a part of the accreting gas can escape, producing an outflow or wind. In the case 
of the ion supported ISAFs, this wind would be thermally driven by the temperature of the 
ions, like an 'evaporation' from the accretion torus. In the case of the radiation supported tori, 
which exist only at a luminosity near the Eddington value, but with much lower temperatures 
than the ion tori, winds driven by radiation pressure could exist. 

The possibility of outflows is enhanced by the viscous energy transport through the disk. 
In the case of thin accretion disks (not quite appropriate in the present case, but sufficient 
to demonstrate the effect), the local rate of gravitational energy release (erg cm^^s^') is 
W = Evrdr{GM/r). The local viscous dissipation rate is {9/4-)vEQ^. As discussed in section 
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|5.4| they are related by 

ev-3[l-(^)i/>, (92) 

where ri is the inner edge of the disk. Part of the gravitational energy released in the inner disk 
is transported outward by the viscous stresses, so that the energy deposited in the gas is up to 
three times larger than expected from a local energy balance argument. The temperatures in 
an ADAF would be correspondingly larger Begelman and Blandford (1999) have appealed 
to this effect to argue that in an ADAF most of the accreting mass of a disk might be expelled 
through a wind, the energy needed for this being supplied by the viscous energy transport 
associated with the small amount of mass that actually accretes. 

Most outflows from disks such as the observed relativistic jets (and the non-relativistic 
ones in protostellar objects) are now believed to have a magnetic origin, requiring the pres- 
ence of a strong, ordered (large scale) magnetic field anchored in the disk. Three- dimensional 
numerical MHD simulations of such processes are becoming increasingly realistic (e.g. Kro- 
lik & Hawley 2010, Moll 2009). The origin of jet-friendly magnetic field configurations, on 
tge other hand, and the reasons for their apparently unpredictable presence are still unknown 
(not all black holes show jets, and the ones that do don't have them all the time). There is, 
however, a degree of correlation of their presence with the hard X-ray states in these objects 
(in the sense that they appear absent in soft states). For a discussion and interpretation of 
these issues see Spruit (2008). 
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